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ABSTRACT 


A pipeline technique by Loomis and Sinha has been applied to the design of recur- 
sive digital filters. Recursive digital filters operating at hitherto impossibly high rates can 
be designed by this technique. 

An alternate technique by R. Gnanasekaran allows high speed implementation using 
the state-space structure directly. High throughput is also achieved by use of pipelined 
multiplv-add modules. The actual hardware complexity will depend upon the number 
of pipeline stages. 

These techniques are used for the design of the IIR notch filter and finally, a com- 


parison of the performance and complexity of these two techniques 1s presented. 
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I. INTRODUCTION 


A. DIRECT FORM STRUCTURE REALIZATION 

A scheme for high speed recursive digital filter realization was suggested by H. H. 
Loomis and B. Sinha [Ref. 1]. The main aspect of the design is in modifying the filter 
difference equation so that the recursive part requires only the Lth past output instead 
of the immediate past output. This removes the tight recursive requirement thus enabling 
pipelining of multipliers and adders in the direct form implementation resulting in high 
throughput. One of the difficulties with this method is that the modification introduces 
extra poles (not necessarily stable) and zeros in the filter transfer function which cancel 
each other in the infinite arithmetic implementation. In practice due to the finite word 
length effect this cancellation may not occur. This may result in an unstable realization 
even if the original filter were a stable one. However, if L is chosen to be large enough, 
then 1t has been proven that the new poles introduced will also be within the unit circle 
thus guaranteeing stability even if the pole-zero cancellation does not occur. This 


method 1s called the direct form structure, or the scalar structure. 


B. STATE SPACE REALIZATION 

Another method by R. Gnanasekaran 1s an alternate technique which relaxes the 
recursive nature of the state equation (state space realization) [Ref. 2]. That 1s, the 
present state vector computation does not require the unmediate past state vector. The 
proposed technique allows high speed implementation using the state-space structure 
directly. High throughput is achieved by use of pipelined muluplv -add modules. The 
actual hardware complexity will depend upon the method of pipelining. though very high 
speed is possible bv incveasing the number of states in the pipelining. This method 1m- 
plements A“ instead of A. the system matrix. The statespace realization will be stable 
if the original filter is stable. Furthermore, the A? filter has improved tolerance to coef- 


ficient error. 


C. OVERVIEW OF THESIS 

Chapter II describes the IIR notch filter. The transfer function and the coefficients 
that are appropriate for the IIR notch filter are developed. 

Chapter III presents the realization of the high speed IIR notch filter by using the 


pipeline technique of Loomis and Sinha. Stability testing, frequency response and unit 


sample response are described. System designs for pipelined multiply-add modules are 
included. 

Chapter IV presents Gnanasekaran’s state space realization which includes a modi- 
fied state space equation. Svstem designs for pipelined multiply-add modules are de- 
scribed. 

Chapter V contains a comparison of the results of the Loomis/Sinha pipeline re- 
cursive design and the Gnanasekaran state space realization design. 

Appendix A contains the program and data for frequency response analysis of the 
IIR notch filters in Chapter III. 

Appendix B contains the program and data for unit sample response analysis of the 
IIR notch filters in Chapter ITT. 


to 


Il. TR NOTCH FILTER DESIGN 
From Rowe, [Ref. 3: p. 72] the IIR notch filter has the transfer function 


byz” + bz + by 


H(z) 
Qyz* + bz = a> 


(2.1) 


Dividing the numerator and denominator of (2.1) bv a, we have 





Do 


b, as a 
a, = Bo Go = A, and -y- = Ay. Therefore the transfer function 1s 
0 0 0 


Let 


Boz" + Bie + Bo 


H(z) = — 
2 + Det Ay 


iaaonm ©.) ang (2.5) we can See that 
1. B, in the numerator, must be equal to B, in the denominator; and, 


Pete ist G, in thie numerator must be equal to the last B, in the numerator. 


Rewriting (2.3) as 


yi B, 
Bz ar +1) 
a) — — (2.4) 
4 ae Byz = Ag) 


the transfer function H(z) is of the form 
Bo{z — 1\)(z — ny) 


eae — 2,) 


(2.5) 


Let the zeros be placed at le=*, and let the poles be placed at A, ex®. Therefore 


the transfer function 1s 


Beales eee") 


OO —————— 
Ce Ome Se) 
H(c) Boz? —let?z - le "2 41] 
D2 Ae, — aes + JA 
+j9 , -j0 
ea € 
| WS = Il) 
H(z) = a 
[z0 — 2 do (Jz + A] 


Bolz? —2 cos Gz +1] 


H(z) = ; ae 
lig a EEG COS Gz + Ao] 
Boz” = 2Bo COs Gz + By 
ies ee 


ae as 2 Apo COS Oz az Ao 
According to (2.3) we can see that 
—2B, cos 0z = —2,/Ay cos 8z 


which means 


For our example of the I1R notch filter we want the zeros to be placed at 


J et/0-4467110339% _ | 4/0 


and we want the poles to be placed at 


+/0.44671103392 a iO 
0.99 ese ae” — ice 


We can see that 
@ = 0.4467110339x 


Ao = 0.9801 


By =/A4n = 9.99. 


(2.6) 


Qa) 


(2.8) 


(2.9) 


(2.10) 


@ip 


(2.12) 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


This realization is shown in Figure 1 on page 6. 


Pym equation (2.1 °) 


4 
Boz” — 2By cos 6z + By 


zt 2,49 Cos Oz + Ap 
Substituting 8, B, and A,. we then have 


0.992* — 2(0.99) cos(0.4467110339z)z + 0.99 


H(z) =— a 
z~ — 2(0.99) cos(0.4467110359z)z + 0.9801 
H-) = 0.9927 — 0.329932 + 0.99 
ee wes — 05 2995- +0.9801 
6 Aa ee ae _-2 
H(27') = 0.99 —0.32993z § +0.99z 


Means | 0.980le 


Because we have the transfer function (Ref. 1: p. 273] 
H =) 4 ic _ 7 a= - 
ie ) —_ a Ae 


from (2.21) and (2.22) we can see that all of coefficients are 


AQ = 0.99 

a = —0.32995 

a. = 0.99 

B, = 0.52995 

po —O9501 . 
= ¥(=) 

oo") 


¥(z) _ 0.99 — 0.3299327' + 0.9927" 
Mel 0 329968 + 0.980ia 


Y(2)[1 — 0.329932! + 0.9801z77] = X(z)[0.99 — 0.3299327! + 0.99277] 


(2.18) 


(2.19) 


—~ 
to 
tJ 
Loy 

wma’ 


(2.26) 


¥(2) — 0.3299327! (=) + 0.980127? F(z) = 0.99.X(z) — 0.3299327 1 X(z) +. 0.99277 A(z) (2.27) 


{fy 


8 = 0.44671103397 





Figure 1. Poles and Zeros Location 


¥(z) = 0.99.X(z) — 0.32993z7 LY(z) + 0.99z-7.X(z) + 0.329932 | ¥(z) 


2228 
— 0.980127 ¥(z) rs 


s(n) = 0.99x(n) — 0.32993.x(n — 1) + 0.99x(n — 2) + 0.329935(n — 1) 


, 
~ 0.9801j(n — 2). ane) 


This is the basic second-order recursive filter (p= 0). The number of p will be ex- 
plained later in Chapter III and the realization is shown in Figure 2 on page 7. 

We can design this IIR notch filter to be a high speed IIR notch filter by using 
pipeline techniques so that it can be operated at hithero impossibly high rates. In 
Chapter III] we will apply Loomis’Sinha method and in Chapter IV Gnanasekaran’s 
state space realization. A comparison of the results of these two techniques will be 


shown in Chapter V. 





= Unit Delay 


= Pipeline Add Unit 


x(n) = 


y(n) = Output 


+> = Pipeline Multiplier 


Input 


Figure 2. | Second Order IIR Notch Filter: p = 0 


Hil, REALIZATION OF HIR NOTCH FILTER BY USING THE 
LOOMIS/SINHA METHOD 


A. PIPELINE RECURSIVE DESIGN 
From Chapter II we have the second order IIR notch filter that is described by the 


transfer function 


(a + we @,z~”) 


H(z7') (3.1) 
(25 =e) 

and the difference equation 

y(t) = O31 — 1) + bay (4 — 2) + agx(n) + ayx(n — 1) + aax(n — 2). (323) 


The computation of 4(77) requires the immediate past output }(m — 1). This filter can 
be developed into high speed IIR notch filter by using the pipeline technique of Loomis 
and Sinha [Ref I]. The main aspect of the design is in modifving this filter difference 
equation so that the recursive part requires only from the Lth past output instead of the 
immediate past output. This means the speed linuts can be overcome by modifving this 
difference equation so that ) (77) depends not on the immediate past output, }(7 — 1), but 
on y(7 — L) for some L > |. This can be accomphshed as follows: 


For L = 1 the difference equation is 
V(t) = by — 1) + byv(m — 2) + ayx (at) + + ayx(n — 1) + @x(n — 2). (3:3) 


We can see that this equation is a 2-order difference equation. 

From [Ref. 1: p. 274] the difference equation has (p + m)"-order where m is the order 
of the difference equation and p is the pipeline delay used to overcome the speed hmit. 
Therefore equation (3.3) has the number of p=0. We can call this equation the second 
order IIR notch filter with p=0. 


Delaying equation (3.3) by 1 yields 
y(t — 1) = byy(n — 2) + byy(n — 3) + agx(n — 1) + a, x(n — 2) + ayx(n — 3). (3.4) 


For L = 2, putting (3.4) into (3.3) results in the difference equation 





3(v) = 6, Lbyy(n — 2) + b(n — 3) + Ayx(i — 1) + ayx(a — 2) + ayx(n — 3)] 


+ byy(nt — 2) + agx(n) + ayx(m — 1) + ayx(n — 2) = 
¥(”) = bey(n — 2) + by byy(n — 3) + agb x(n — 1) + aybyx(n — 2) + anh x(n — 3) 3.6) 

+ byy(n — 2) + agx(n) + ayx(n — 1) + ayx(n — 2) - 
y(n) = (b? + by)y(n — 2) + b)byy(1 — 3) + agx(n) + (agb; + a,)x(n — 1) 3.7) 


+ (a,6, + ay)x(n — 2) + ayb, x(n — 3). 


We see that the computation of »() requires only the second past output y(n — 2) 
instead of immediate past output 3(#—1), and we can use this technique for the 
higher-order of L. 

We see that equation (3.7) is a 3-order difference equation and has p=1. Thus we 
can call this filter the augmented IIR notch filter for p= 1. 


The z transform of equation (3.7) is 


eben + bbz VY ta X¥ + (qb, +a)2 X 


3.8 
(a,b, + a)z -X + a5b,z °X (3.8) 
. 1 — (bj + by)2* - bbz” 
=I ag) a 
= Ag taz +az “\Il+ 0,2 
PP ida Fae NI + ye) a 


ae lab) 


That means the transfer function of the augmented IIR notch filter, p= 1, 1s equal 
to the transfer function of the second order IIR notch filter, p= 0, multiplied and divided 
by the same polynonuals. 

We can design the higher-order difference equation by starting from a second-order 
difference equation. That means we can design the higher order, p-augmented IIR notch 
filter by starting from our second order IIR notch filter, p=0, as follows: 

From Chapter ii the second order IIR notch filter transfer function is the original 


transfer function 


_ (a+ Ge Dee 


H(z7') . 
(]— bz — 427’) 


(3.11) 


Multiply the numerator and denominator of the original transfer function by the 
same polynomial equation. This can be thought of as equation (3.12) As developed in 
(Ref. 1: p. 281]. the transfer function of the p-augmented filter is as shown rl equation 
(3.13) We know from other considerations that the transfer function of the higher order 
realization must in fact be the same as the original and therefore must have the same 
poles and zeroes in the z - plane. 

The coefficients of the a polynomial depend on the original 5, and 6, and p is de- 


termined by the number of stages in the pipeline multiply and add units. 


(a) + az =F doz )(1 +f o2—! Seca a ag em) 


H'(z7') = (3.12) 
(l— bz! — byz~*)(1 ~ o2—" +o + alee) 
We ay teh aye i 7 
er ear ia _ 
prit p+2° 


The values of 5%, a. and a? are given in Tables 1, 2 and Ret Wp 275-2 


Table 1. s for oes = 


bih, + 2b,b? 





es : | 
4 bs + Abb, + 3b,b3 bib, + 3b2b3 + b3 
5 bo + Sb4b, + 6b2b? + B3 bib, + 4b3b? + 30,63 


10 


Table 2. « 'sforp = 0.1.. 


| 
EE ee es 
DE ee ee 
2 it [he |e+e fo fo 
a a a 
Ee ee 
Ee en ee ee 











Table 3. a’sforp = O.1.. 


a 
(te a 


+ ab} + 26,6,) 
ae Ae ea 























nies 20,03) 
+a(bj+ 5) 
+ aD, 










+ a,(b} + 2b,/,) 


B. SECOND ORDER IIR NOTCH FILTER: P = 0 


From equation (3.13) in section A, we let p = 0. This equation becomes: 


—| —2 
(a + Aye -+- M52 ) 





H'(27') = — = (elt) 
(l—6,z  —65z *) 
a Eta) 

a 1 = \ a 
ee (z ) Nz (3 15) 

: _ =" 

a - (dy + a2 a5 a) (3.16) 
A(z) (1 — bz! — bz”) 
Manipulating (3.16) will result in the difference equation, (3.20). 
¥iz)E1 — b27! = by277] = A(2)Lag + yz) +277] (3.17) 


J) 


¥(z) — byz7) ¥(z) — byz 7? ¥(z) = ag X(z) + ayz7'X(z) + 27° X(z) (3.18) 
¥(z) = agX(z) + ayz7 UM) + ayz77N(2) + Bz! ¥(z) + 527° ¥(z) (3.19) 
y(n) = Apx(m) + ayx(m — 1) + anx(n — 2) + by(n — 1) + dyn — 2). (BPA) 


This realization is shown in Figure 3 on page 13. 


From chapter II we know that for our example notch filter, that the coefficient 


values are 

aA = 0.99 

a, = —9.32993 

ay = 0.99 (3.218 
b, = 0.32993 

6, = —0.9801. 


The transfer function 1s 


= =O 

_ Aye 2 ee ee) 
H'(z Jepccnntlais (3.22) 

(l-—b,z  —h,: 7) 

The difference equation 1s 
¥(2) = dgx(4t) + Qax(u — 1) + agx(a — 2) + by (mn — 1) + b33(7 — 2). (3225) 


C. AUGMENTED IIR NOTCH FILTER: P = ! 


From equations (3.12) and (3.13) in section A, we let p = 1. Those equations be- 





COME, 
H'(27") = (aq tqz 7 taz V+ ay") (3.24) 
CUS tye =e ON shee |) 
(1). () -t | (1) -2 , 1) ,-3 
e Ay ay 2) a ae 
came 5) wa = OF Zz 
re Jha) : 
H'(z7') = To (3.26) 





Figure 3. Second Order IIR Notch Filter: p = 0 





Peea. + a,'z + ae oe afz-4 
TO) — 3. = (3.27) 
— b, ~ — b; a 
Manipulating (3.27) will result in the difference equation, (3.31). 
Mee, 2s — b, 27 = X(2)lay + a,’z + as 2° + ay ’2 7) (3.28) 
H(z) — 662-2 F(z) — OM 2-2 (2) = al) + aMM27 Lz) + al2-2X(2) oe 
os 
+ as2~7 Xz) 
¥(z) = af) X(2) + afP27 UN (2) + af)27*X(2) + af 277 Xz) + 8372" Y(z) aa 
4 Le Y(z) | 
y(n) = a x(n) + al x(n — 1) + a x(n — 2) + a x(n — 3) + BSP (0 — 2) ay 
=F a y(n — 3). 8.5) 


This realization is shown in Figure 4 on page 14. 
From Tables 1, 2 and 3 [Ref. 1: pp. 278-279], we have 
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Figure 4. The Augmented IIR Notch Filter: p = 1 


an! = | 

a!) = b = 0.32993 

ay? = ay = 9.99 

at” = ayb, + a, = —0.0032993 

as = a,b, + a, = 0.8811462 

a3’ = ayb, = 0.3266307 

bs”) = bf + by = —0.8712462 

b§? = bb, = —0.32336439 . 
The transfer function is 

ae = a)z—! ar a§)2~? a are 


es ae 
i. ]—- ae — E),- 


The difference equation is 


(3.32) 


(3.33) 


(2) = ah x(n) + a} x(n — 1) + a8 x(n = 2) + a x(n — 3) + BO ¥(n = 2) 


ee viru 3). 


D. AUGMENTED IIR NOTCH FILTER: P = 


From equations (3.12) and (3.13) in section A, we let p = 


come: 
H'(c7') = (a + oe + ay2z~°\(1 -- of) 2! “te a2?) 
> an —} — Oe = 
(l—bz —boz “\lt+ay’z  +a5%2z *) 
Hoe") — ag + qy)z7} + az? + a@,~ + a. 
= a ee 
3 z 4 Ae 
¥(z) 
H(z‘) =—— 
A(z) 
bz} a? a te 2). Ti = ee = a. rE gQ@,— 
ie (22,2 Spe, 


Manipulating (3.38) will result in the difference equation, (3.42). 


Y(z)(1 — bz? — E244 = X(z)[a” _ aze + a. + az? 


¥(z) — bY2-? ¥(z) — oz Wz) = ak ve )tatz es roy ela 
2 az 3¥(2) + ay iz ~* X(z) 


Y(z) = aX (2) + az X(2) + a2 ON (2) + a2 OA X(z) + a2 F(z) 


+ BF ¥(z) + 62 ¥(z) 


in) = a x(n) + aM x(n — 1) + a x(n — 2) + a x(n — 3) + a? x(n — 4) 


SD. in — 3) + OS y(n — 4). 


This realization is shown in Figure 5 on page 16. 


From Tables 1, 2 and 3 [Ref. 1: pp.278-279], we have 


(3.34) 


2. Those equations be- 


(3.39) 


(3.40) 


(3.41) 


(3.42) 





Figure 5. The Augmented IIR Notch Filter: p = 2 


om = 1 

og = b, = 0.32993 
2 = b; + by = —0.8712462 
ay? = = a, = 0.99 

a\” = agb, + a, = —0.0032993 


(3.43) 
= ay(b? + by) + a,b, + a = 0.01861246 
a = a,(b; + by) + anb, = 0.61408096 
ay” = a,(b; + by) = —0.86253373 
bi) = b* + 2b,b, = —0.61081465 
bY = b?b, + bs =0.8539084 . 
The transfer function is 
Q) 4, Qo 14 42)-2 4 .2).-3 4 2),-4 
zt ap a OF or Q3 “a 
oS a (3.44) 
i =a 
The difference equation 1s 
ie?) — a? x( 1) + a x(n —1)+ ay x(n — 2)+ a x(n — 3) (3.45) 
45 


a x(n — 4) + BS y(n — 3) + bP y(n — 4). 


FE. AUGMENTED ITR NOTCH FILTER: P = 3 


From equations (3.12) and (3.13) in section A, we let p = 3. Those equations be- 





come: 
mil =o 3 yess 
Hh = (a+ az taz “\l+ a pees es, =e ats es 3.46) 
eo) = =i =) , 3) 3),-2 G3 spake 
(1 =e 4 Soa )Q + a, z ee + 03 ra ) 
a a ae az + a2-? - al? ee eae ey oe ae2- : 
H(z )= Cee oe (3.47) 
l = b ZL = bs Zz 
Pe ¥(z) 
H'(z7') = Wo (3.48) 
¥(z) _ &. 1 azo! — a z- ee a), + a2" + az (3.49) 
os os 
A(z oe = ee 


7 


Manipulating (3.49) will result in the difference equation, (3.53). 
¥(2)[1 — 63? — 62>] = 
2) 


MM Y(z)Cay” a fo z+ aS2~* oF az” gs af a 


¥(z) — B27 ¥(z) — BP’z-? Vz) = af? X(z) + afz7'X(z) + a2? X(z) 
+ az X(z) + az X(2) (3.51) 
— a X(z) 


¥(z) = aS X(z) + a2 X(z) + a)z77.X(z) + az X(z) + af’ zF7X(z) 


z= 2 (3:52) 
+h ae a A(z) + ba ¥(2) + bP) > ¥(z) 
y(n) = as (1) + a? x(n = 1) + asx (n — 2) + a? x(n = SP + a? x(n _— 4) G3 53) 
Joe 
+ a(n — 5)-+ by Gee ge): 
This realization is shown in Figure 6 on page 19. 
From Tables 1, 2 and 3 |Rel_ 1: pps 2732279) citar 
a?) = | 
ag) = b, = 0.32993 
a8) = b? + by = —0.8712462 
a) = b} + 2b, b, = —0.61081465 
as = ay = 0.99 
a?) = ayb, + a, = —0.0032993 | 
(3.54) 


a) = ap(b? + by) + a,b, + a, = 0.01861246 
ay) = ao(b} + 2b, by) + a,(b} + by) + a,b, = 0.00937445 
a? = (bi + 2b) by) + ay(b? + by) = —0.66 100766 
ae? = a(6j + 2b,b,) = —0.6047065 
po) — by + 3b7b, + by = 0.65238232 
bY) = bi} b, + 2b,b3 = 0.59865944 . 


The transfer function is 


an) gee ay 2a de ia ae 


77) Bessa aii a are a 
H'(27') = B= 53 | (3.55) 


The difference equation is 
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Figure 6. The Augmented IIR Notch Filter: p = 3 


y(n) = as x(n) oP ee —I)+ aS x(n — 2) 
+ as x(n — 3) + at x(n — 4) 


+ al x(n —5)+ b? y(n — 4) 


aL bP y(n — 5). 


(3.56) 


F. STABILITY TESTING 

From section A if we start with a stable filter, we would like to be assured that the 
poles of the augmented-order filter all be outside the unit circle in the z~' plane. We de- 
sired this because, even though the added poles are cancelled by the zeroes of 
(1 + ¢@z-! + +» + o6)2-*), realization imperfections will prevent exact cancellation and the 
augmented filter would be unstable. In order to examine the stability question, we will 
use Jury’s stability test as applied to successively higher-degree denominator polynomials 
in Z as p increases [Ref. I: p. 281]. The results from this paper are summarized below: 

1. Condition for Stability of the Original IIR Notch Filter: p = 0 

The original IIR filter is described by equation (3.14) and to test for stability, 


we analvze the following characteristic polynomial resulting from it: 
Z: 
F(z) =2° -b,z-b,=0 (32090) 


The conditions for stabilitv are: 


a ee) 
1+5,-b,>0..9 (3.58) 
[2 i) 


These conditions are shown graphically in Figure 7 on page 21. 
Circled numbers in the figure refer to the conditions of the same number in equation 
(S508) 

Checking for stability of our example IIR notch filter; we have from © of (3.58) 


1—b,-—b,>0 
1 — (0.32993) — ( —0.9801) > 0 
1.65017 > 0 


from @ of (3.58) 


NM 





Figure 7. Condition for Stability of the Original Filter [Ref. 1: p. 283] 


110152993) — ( —0.9801) > 0 
2.31003 > 0 


from ® of (3.58) 


Heal | 


10.9801] <1 
0.9801 <] 


Thus the conditions are all satisfied and the original filter is stable. 
2. Condition for Stability of the Augmented IIR Notch Filter: p= 1 
The augmented IIR notch filter, p= 1, is described by equation (3.33) and to test 


for stability, we analyze the following characteristic polynomial resulting from it: 
F(z)=z+6,=0 se) 


The conditions for stability are: 


21 


1+5,>0..0 


i-b,>0.,. 9 Oe 


These conditions are shown graphically in Figure 8 on page 23. 
Circled numbers in the figure refer to the conditions of the same number in equation 
(3.60) 

Checking for stability of our example IIR notch filter; we have from © of (3.60) 


1+6,>0 


1+ 0.32993 > 0 
eo 2295 20 


from © of (3.60) 


0.67007 > 0 


Thus the conditions are all satisfied and the augmented IIR notch filter, p = 1,15 stable. 
3. Condition for Stability of the Augmented IIR Notch Filter: p = 2 
The augmented IIR notch filter, p= 2, 1s described by equation (3.44) and to test 


for stability, we analyze the following characteristic polynomial resulting from it: 
F(z) =2° + bz + (b; + bs) (3.61) 


The conditions for stability are: 


1+b, +(b; +5,)>0... 0 
loki Eee. & (3.62) 
aA Sb eis. 8 


These conditions are shown graphically in Figure 9 on page 24. 
Circled numbers in the figure refer to the conditions of the same number in equation 
(3.62) 

Checking for stability of our example IIR notch filter; we have from © of (3.62) 


1b by SC 
1 + (0.32993) + [(0.32993)? + ( —0.9801)] > 0 
0.4586838 > 0 





Figure 8. Condition for Stability of the Augmented Filter: p = 1 [Ref. 1: p. 284] 


from © of (3.62) 


(0152993) + [(0.32993)° + ( -0.9801)] < 0 
—0.2011762 <0 


therefore the filter is unstable, from © of (3.62) 


ee 


| (0.32993)? + (—0.9801)| <1 
0.8712462 <1 


Note that @ is violated and here the realization will be unstable. 
4. Condition for Stability of the Augmented IIR Notch Filter: p = 3 
The augmented IIR notch filter, p= 3, is described by equation (3.55) and to test 


for stability, we analyze the following characteristic polynomial resulting from it: 


F(z) =z? + bz? + (b> + by)z + (b} + 26,62) (3.63) 


23 





Figure 9. Condition for Stability of the Augmented Filter: p = 2 [Ref. 1: p. 285] 


The conditions for stability are: 


ee om 2b (| So Ub = O., © 
1 —b, + (b; + b,) — (6° + 26,5,) >0... @ 
|b? + 2b,5,| <1... @ 

(62 + 2b,b,)? — 1] > | b)(b? + 2b.) — (62 +5)... © 


(3.64) 


These conditions are shown graphically in Figure 10 on page 25. 
Circled numbers in the figure refer to the conditions of the same number in equation 
(3.64) 

Checking for stability of our example IIR notch filter; we have from © of (3.64) 


ak by bee Se Epa Gy 0) eo 


1 + 0.32993 + [(0.32993)" + ( —0.9801)] + [(0.32993)? + 2(0.32993)( —0.9801)] > 0 
~0.15213085 <0 
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Figure 10. Condition for Stability of the Augmented Filter: p = 3 [Ref. 1: p. 286] 
therefore the filter is unstable, from @ of (3.64) 


L2G SSS DS 
1 — 0.32993 + [(0.32993)” + ( —0.9801)] — [(0.32993)° + 2(0.32993)( —0.9801)] > 0 
0.40963846 > 0 


from ® of (3.64) 


|b? + 2b,b,1 <1 
| (0.32993)? + 2(0.32993)( —0.9801)| <1 
0.61081465 <1 


from ® of (3.64) 
on 3 2 ’ 3 2 
1 (by + 2b,by)° — 1] > 164(b; + 2,62) — (b; + 62) 


| (0.32993)? + 2(0.32993)( —0.9801)}? — 1] < 


| (0.32993){(0.32993)° + 2/0.32993)( —0.9801)} — {(0.32993)* + ( —0.9801)} | 
0.62690546 < 0.66972012 


ae 


therefore the filter is unstable. 


Note that © and © are violated and therefore the realization will be unstable. 


G. FREQUENCY RESPONSE 

From the transfer functions (3.22), (3.33), (3.44), and (3.55), we use the program in 
Appendix A to calculate the magnitude and phase of the frequency response of the se- 
veral realizations. We plot the magnitude of these responses for the original filter and for 
the pipeline augmented that all of the frequency responses are the same from those fig- 
ures. 

1. Second Order IIR Notch Filter: p = 0 

The magnitude response is shown in Figure 11 on page 27. 
The Augmented IIR Notch Filter: p = 1 
The magnitude response is shown in Figure 12 on page 28. 
3. The Augmented IIR Notch Filter: p = 2 

The magnitude response is shown in Figure 13 on page 29. 
4. The Augmented IIR Notch Filter: p = 3 


The magnitude response is shown in Figure 14 on page 30. 


to 


We know from section F that the example augmented IIR notch filter for p = 
2, and 3 are unstable but we cannot see that from these magnitude responses. We have 
to find another way to show which one is stable or unstable. This will be shown in sec- 


oral JL. 


H. UNIT SAMPLE RESPONSE 

From the difference equations (3.23), (3.34), (3.45), and (3.56), we use the program 
in Appendix B to calculate the output for a unit sample input to each of the pipeline 
realizations. We plot the outputs in Figures 15, 16, 17 and 18. We can see that the out- 
puts for Figure 15 (p=0) and Figure 16 (p= 1) are stable, however the outputs for Fig- 
ure 17 (p= 2) and Figure 18 (p= 3) are unstable as predicted by the stability testing in 
Section i. 

1. Second Order IIR Notch Filter: p = 0 
The output is shown in Figure 15 on page 31. 
The Augmented IIR Notch Filter: p = | 
The output is shown in Figure 16 on page 32. 
3. The Augmented IIR Notch Filter: p = 2 


The output is shown in Figure 17 on page 33. 
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Magnitude Response of Second Order IIR Notch Filter 


Figure 1. 
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Figure 12. 
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MAGNITUDE RESPONSE 
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Figure 13. 


a.0f 


x 10xe 8 


gw 
tog Oe 
| 
Cw oe | 
aw a | 
al 
— 
— 
| 
) Ee | 
| Eo | 
| ne | 
| an. | 
= 
5 mJ or aor fn poe 
on as =z ow | aa =e 
| aa? | | on | cw | n> | Cc 


MAGNITUDE 


Magnitude Response of the Augmented IIR Notch Filter: p 
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Magnitude Response of the Augmented IIR Notch Filter: p 


Figure 14. 
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Output Sequence of Second Order IIR Notch Filter: p 


Figure 15. 
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4. The Augmented IIR Notch Filter: p = 3 


The output is shown in Figure 18 on page 35. 


I. PIPELINE IMPLEMENTATION 

All of the IIR filters in sections B, C, D, and E can be developed into high-speed IIR 
filters by applying pipeline techniques. We use Ka-stage pipeline add units and Km-stage 
pipeline multipliers to increase the speed of the filters. The Km stages are involved in the 
coefficient multiplies and the Ka stages are involved in the adds. The increase in order 
of the filteris P= 2Ka+ Km— 2 |Ref. I: py 280). 

1. Second Order IIR Notch Filter: p = 0 

To achieve a p=0 realization, we cannot use pipeline modules with longer de- 


lays than Aa= 1 and Am = 0. 


P=2Ka+ Km—2 
=2(1)+0-2 
==) 


We have a realization of the original IIR notch filter but with an additional 3 


delays. The equation is 
z°Y= Baal cae Syn ee eC tiie i ble Ya (3.65) 


This realization is shown in Figure 19 on page 36. 
2. The Augmented IIR Notch Filter: p = 1 
To achieve a p=1 realization, we cannot use pipeline modules with longer de- 
lavs than Aa@=1 and Am = 1. 
P=2Ka+Am-2 
=2(1)+1-2 
=] 


We have a reahzation of the augmented IIR notch filter, p = 1, but with an 


additional 4 delavs. The equation is 
z~ Yaz [aX + ajz X+4 az X4 az X + boz “Y + 6,2 “Y|. (3.66) 


This realization is shown in Figure 20 on page 37. 
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Figure 19. Second Order IIR Notch Filter: p = 0 with Ka=1, Km=0 


3. The Augmented IIR Notch Filter: p = 2 
To achieve a p=2 realization, we cannot use pipeline modules with longer de- 


lavs than Ka= 1 and Am = 2. 


We have a realization of the augmented IIR notch filter, p = 2, but with an 


additional 5 delays. The equation is 
2 a PN one == aM Ss az X + ace + ae b,z~° y + biz val (3.67) 


This realization is shown in Figure 21 on page 38. 
4. The Augmented IIR Notch Filter: p = 3 
To achieve a p=3 realization, we cannot use pipeline modules with longer de- 


lavs than Ka=1] and Km=3. 
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Figure 20. The Augmented IIR Notch Filter: p = 1 with Ka=1, hm=1 
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2 with Ka= 1, Km=2 


nted IIR Notch Filter: p = 


The Augine 


ae 


Figure 


We have a realization of the augmented IIR notch filter, p = 3, but with an 


additional 6 delavs. The equation is 


x - X+Qaz_ 'X ee X PAO 
ees |° | (3.68) 


mie Cie Seay Nase bag Se Yeeer oY, 


This realization is shown in Figure 22 on page 40. 


Table 4. Saanioamem PIPELINE RECURSIVE DESIGN 





J. SUMMAIARY 

This chapter investigated a method for applving pipeline techniques to the design 
of a high-speed IIR no’ch filter. Using the Loomis Sinha method, IIR notch filters op- 
erating at rates hitherto impossible can be designed. The general structure of the IIR 
notch filter and the ™2thod of calculating the multiplier coefficients have been presented. 
The stability of the resulting realization has been investigated and a technique for as- 
certaining the stability of the realization has been presented. The delay and complexity 


of each realization are summarized in Table 4. 





Figure 22. The Augmented IJR Notch Filter: p = 3 with Ka=1, Km=3 
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IV. GNANASERKARAN’S STATE SPACE REALIZATIONS 


IIR filter realizations are possible and most of them can be presented within the 
general time domain representation of the system under consideration. In continuous 
time It Was convenient to represent a linear constant coefficient differential equation by 
a system of first-order differential equations for analysis purposes. It is, therefore, not 
surprising that it is useful to represent a linear constant coefficient difference equation 
by a system of first-order linear constant coefficient difference equations. Such a for- 
mulation results in a special form of representation of the difference equation called a 


state variable representation. 


A. STATE SPACE REPRESENTATION OF A NTH ORDER RECURSIVE 
DIGITAL FILTER 


We are used to working with a difference equation of the form 
y(n) = b(n — WI) + + bay (n — VN) + agx(t) + aqax(n -— I 4 teapx(n— MY (4.1) 


or in the notation 


NV M 
v= > bn — b+ Y ayx(n — K) (4.2) 
c=) k=0 


that is called the Nth-order form. An alternative is to describe an LTI system with N 
first-order difference equations and an output equation. This alternative form for de- 
scribing a linear time-invariant system is known as the state-space description and can 


be Written as a matrix equation: 

v(1) = Av(n — 1) + Bx(x) (state equation) (4.3) 
with the system output given by 

(4) = C’v(n) + dx(n) (output equation). (4.4) 


In Eqs. (4.3) and (4.4), v(7) is a vector with N rows and | column (an .V x I vector or 
column matrix) of newly defined variables called the states, x(7) is the system input, and 


y(n) its output. In the matrix state equation of Eq. (4.3) the matrix A must be an 


4] 


Nx N matrix and the B matrix an Nx 1 vector. The scalar output y(z) 1s defined by Eq. 
(4.4) where C is a Nx 1 column matrix and d is simply a real number that multiplies (or 
scales) the input x(7). The elements in A, B, C, and d are constants that depend only 


upon the system parameters. 


B. WRITING THE STATE SPACE EQUATIONS 
To illustrate this approach to modeling of LTI systems we return to a version of the 


Nth-order difference equation described as a second order IIR filter, namely 
Wert) = dyn — 1) + boy = 2) again) + Gxt — 1) aax(n = 2) (4.5) 
Y= be ¥ oy) Vi ee eee (4.6) 


We define the state varible as: V,, 14. Then, we can establish the following re- 


lationships: 


MH Vs Digs 0 bz" + Gh eee er anz* X 
Voez Vi =27'Y (4.7) 


1), 2 
Zz VM, = Z Y. 
Therefore the state equation ts: 


Vi = bz) Foe been aXe ene a wk 


ei, (4.8) 

We express it in the matrix state equation: 
v(72) = Av(n — 1) + Ba(n) (4.9) 
V=Az V+BX (4.10) 


V b, b V Gy Etna 
] = ] Z o- ] I 0 1¢ 2¢ Oe (4.11) 
vi ti of |y, 0 


V 
v=] (4.12) 
V 2. 


That means: 


The output equation will be: 


y(n) = Cv(n) + dx(n) 
Vy 
ashe, eal 7 + dX 
L 


Y=c,V, +c), 4+ dX. 


C. MODIFIED STATE SPACE EQUATION 


From the state equation 


v(11) = Av(n — 1) + Bx(n), 


(4.13) 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


the computation of v(n) requires the immediate past state vector v(n-1). This recursive 


nature severely limits the speed that can be achieved by direct implementation of the 


State equation. The speed limits can be overcome by modifving this state equation so 


that v(n) depends not on the immediate past vector, w(#n—1), but on s(n — L) for some 


L>1. This can be accomplished as follows: 


mon?) — | the State equation 1s 
v(41) = Av(n — 1) + Bx{(n) 
V=Az V+BX. 
Delaving equation (4.19) by 1, 2, and 3 yields 
vin — 1) = Av(n — 2) + Bx(n — 1) 
v(1 — 2) = Av(n — 3) + Bx(n — 2) 
v(n — 3) = Av(n — 4) + Bx( — 3). 


For L = 2 putting (4.21) into (4.22), the state equation, results 1s 
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v(77) = ALAV(n — 2) + Bx(z — 1) ] + Bx(n) 
v(1) = A’v(n — 2) + ABx(n — 1) + Bx(v) 
V = A‘z°V + ABz X+ BY. 
For L =3 putting (4.22) into (4.25), the state equation, results is 
v(n) = A’ Av(n — 3) + Bx(n — 2)] + ABx(n — 1) + Bx(x1) 
v(n) = A°v(n — 3) + A°Bx(n — 2) + ABx(n — 1) + Bx(x) 
V = A’z °V + A’Bz 7X + ABz |X + BY. 
For L = 4 putting (4.23) into (4.28), the state equation, results is 
v(n) = AOL Av(n — 4) + Bu(t — 3)]] + A°Bx(n — 2) + ABx(n — 1) + Bx(n) 
v(1) = A’v(2 — 4) + A°Bx(2 — 3) + A’Bx(n — 2) + ABx(n — 1) + Bx(12) 
V=A'z2_V + A°Bz ?X + A°Bz “NX + ABz X+4+ BX. 


Equation (4.32) can be written in the following form for anv L> 1: 


L—] 


V =  a'Be“s 
= 
— 


(4.24) 
(4.25) 


(4.26) 


(4.30) 
(4am 


(4.32) 


One of the benefits of the state space realization is that if the original filter (L= 1) 


is stable, all higher order realizations are also stable, and furthermore, A* has improved 


error properties over A. This advantage of the state space method 1s discussed in the 


paper by Lu, Lee, and Messerschnutt [Ref. 4]. 


D. SECOND ORDER IIR FILTER STATE SPACE REALIZATIONS 


From the state space equation of the second order IIR filter in section B and the 


modified state space equation in section C, we can develop techniques for high speed 


implementation of second order IIR filters for several values of L. This section shows 


how to apply these concepts for L = 1, 2, 3, and 4: 
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Beginning with equation (4.20), 
V=Az V+BY (4.34) 
and substituting for V, A, and B from section B yields: 

My = b; by = yy 

V; 1 0 V; 
4,35 
* +az + ale aoe) 

=F 0 De 


Beginning with the output equation, equation (4.17), and referring to Figure 23, 
we note that the output of the multiplier dis dz-2¥. This results is a revised output 
equation of 


eee ase \. (4.36) 
Substituting the state equation into the output equation, results in 


She ane 
SO hyo ue lh 
+ age, Xx (4.37) 
rater 


+ (a5c, + dz X 


Ve the 2 aye 
a byez Y 
+ ayc,Xx (4.38) 
ae a,cjz X 
ap 1 @he, ap d)z~°X. 


To demonstrate equality with the Y of the original IIR filter as expressed in equation 


(3.23) for p = 0, let c, = 1, c,=0, and d=0. Then the output equation is: 
aoa Vt baz ok + oN az x Payee. (4.39) 


By substituting: 
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D, = b1, 


b, = by 
dy = @)) 
a, = a1 
9 Ss 


We can Write this equation as: 


Y= bz Y+ biz Y+ ay4X+ mae ee ayqz°X. 


This realization is shown in Figure 23 on page 47. 


2 L=2 


Beginning with the state equation (4.26), 


2 


V =A’z “V+ ABz 1X + BY 


and substituting for V, A, and B from section B yields: 


| 


I 


i 
2 


V, 
eo 


Hl 


b, by |) by by _—2 My 
Sy cia oe, 


| =, 
7 b, by Ay == Qy2 == QZ me 
1 90 0 


a ara aaz~* : 
4° ; xX 


Me] forth, bb] Lf M, 
V, 7 b, b, a. V, 


n ee are a,b,z7' a ayb,z~* 


= Il ao 
AgtaZ +z 


—] —2 
A,+aZ +4@Zz 
[° vita |r 


by + b, bb, | _»| VY 
b, b, ; Vy 


‘| 


ay + (ayb, + a,)z7' + (a,5, + ay)2~° + ayb,2z~° 


=i aes) = 
Gel 2 ih 2 
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3 


fos 


fx 


(4.40) 


(4.41) 


(4.42) 


(4.43) 


(4.44) 


(4.45) 





L=1 


igure 23. Second Order IIR Filter: 


Beginning with the output equation, equation (4.17), and referring to Figure 24, 
we note that the output of the multiplier dis dz-W. This results is a revised output 


equation of 
Y=.) Peo aoe (4.46) 
Substituting the state equation into the output equation, results in 


¥=[(bj + be, + ye |zV, 
MECH ee cb Sele os 
+ ac, X 
+ [(aob, + a)e, + aye, ]z  X 
+ (a,b,c, + ac, + @ycy)z XX 
+ (ab\c, + ac, + dz? X 


(4.47) 


Y=[ (oj Hb) He fz -¥ 
Be We) oe alee le 
+ ayc).4 
i; | (Gaovee Galen anesll am 
+ (a,b,c, + ajc, + a\cy)z7"X 
++ (Ga0Ney One eet 


(4.48) 


To demonstrate equalitv with the Y of the augmented IIR filter as expressed in equation 


(3.34) for p = 1, letc, =1,c,=0, andd=0. Then the outputequation 1s: 


Y= (b; +b) Y-Pbypj2” Vp ap eenias eter 12 Mee 


- er. (4.49) 
+ (a,b, + a,)z oye A,0)\zZ ae 
By substituting: 
2 
by +b, =b,, 
bb, = 01» 
om, (4.50) 


ajb, + a, =a) 
a,b, + a = a)3 
Ab, = ay4 


Wwe can write this equation as: 
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eee em ian X 


z 7 (4.51) 
ane as ere “OK 


This is the same as the equation (3.34) for p= 1 found on page 15. This realization is 
shown in Figure 24 on page 50. 
a L= 3 


Beginning with the state equation (4.29), 
V = A°’z OV + A’Bz 7X + ABz  X¥+BYX (4.52) 


and substituting for V, A, and B from section B vields: 


Vy b; +b, 66, )) 6, | 3) i 
= Z 


n be +6, 6,6) \| ag t+ az + ye -2y 
eae 
b, b, 0 


(4.53) 
b, 6, \ajtajz0 + aoz~° = 
+ ane 
et 0 
Glenn. aoe 
. gta ta, h 
0 
V, by Coole) bale Ga || al i 
= Z 
le bg aby bby V, 
oe, Eee ee by) Ua: ee : 
dnb, + a,b,2- | + ab,z te 
0) Boa | pes (4.54) 


Ay + vd _ AZ 


ob x 


in +a ee ee a,b)z— : -ly 
=> Fa 


‘ae ees a 
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a [4 (h? + 2b,)  b,(b? + by) 
(|e a Dy bb, : 
ay (anoeranee —- Wanda c,beemaio pay \ze (4.55) 


ae oe (a,b + a,b, + felon + (a,b; + @yb.)z~" Xe 
az + (gb, + a,)z~* + (a,6, + a))2~° ae 


Beginning with the output equation, equation (4.17), and referring to Figure 25, 


we note that output of the multiplier d is dz~+X. This result is a revised output equation 
of 


Fac,V,+0V,+4z X. (4.56) 


Substituting the state equation into the output equation, results in 


= [ (db; + 2b,,)c, + (Of + by)e, Jz; 
+ | (6;8, + b3)e, + bybyc, |2 7 V 


+ Goce. 

"eae ney Coen nen 

+ [ (aby + agby + a,b, + aghey + (agb, + a)cq |2~ A 
+ [ (a,b5 + a,b, + ab,)c, + (a,b, + a Joya 

+ { (a,b} + ayby)c, + aybycy + dz 8 


y = [ (o; 4= 2b, bo)e, a (b: == Delay eae I 


+ [ (bib) + b3)cy + bybyeq |z~} 
+ agc, Xx 


+ [(a)6, + a,)c, + ereulz ae 

+r [ees + dajb, + a,b, + ay)c, + (ab, + aenianex 
+ [ (a,b, + a,b) + ayb,)c, + (a,b, + eon 
+ [ (ab? + ayby)c, + abycy +d ]2~ 


To demonstrate equality with the Y of the augmented IIR filter as expressed in equation 


(3.45) for p = 2, let c, =1, c. =0. and d=0. Then the output equation is: 


5S} 


Y =(b° + 26,552” VY Fb byeeo. 2 9 1 an een eee 
+ [ ao(bp + by) + a,b, + a Jz 7X + [a,(; + by) + agb, Jz X (4.59) 
a; sae 


By substituting: 


b} + 2b, b> = br, 
bib, + by = by 
mos 

A; + a = ayy (4.60) 
ag(by + by) + ab, + a = 443 
ay(b} + b) + aby = ay4 
ay(by + by) = ays 


we can Write this equation as: 


Va) 25 ee em aarp care eae 


. 2 Bae (4.61) 
+ Ay3Z ay + Qy4Z BCS ape: ne 


This is the same as the equation (3.45) for p=2 found on page 17. This realization is 
shown in Fisure 25 onpaeere: 
4. L=4 
Beginning with equation (4.32), 


V=A°7—-V + A Bee YA Dee As ee (4.62) 


and substituting for V, A, and B from: section B yields: 


in 
tJ 






“G 





er 
ofesisiaraial 





fr 


Figure 25. The Augmented IIR Filter; L = 3 


My bi(bj + 2b) by(by + by) |, by] al Y, 
V, De ip b,b, lL cog elas 
b (64 +2b,) b,(bf + d)) [Pag tayz) + az7° | _, , 
7A 
nn ok Dh bb, 0 


1 -2 
|? +b, b, al + az + az foe (4.63) 
1 


+ 


b a, +az )+az~ = 
is : 2 _ 2 fk ly 


=? 
An + az M4 az 
_ [ater +a |r 


V ke 3bib, + bs b,bo(bt + 2b,) | _4| V, 
z 
ie b,(b? + 26,) Be(oqer b.) V, 
dob, (b; + 2b,) + a,b\(b? + 2b,)z7! + ayb,(b? + 2b,)= ‘ “sy 
Fs 


-- 
ag(by + 2) + a,(b} to a ‘+ a(by + by)z- 


9 (4.64) 
AP, SF a) Dee at A,b,z 


a 


Ayb, + ab a Bc a sly 
ag + QZ wets QZ 


|" Bley seers fn + Pade ‘+ an(b} as ees 
L 
| 


& 
A + e. Sas i 


V, by - 3b,b, + Ds Dee) ee Le 
im 7 ew + 26) b(bj + by) } ¥ 
dy + (ab, + a)2~' + (agby + Agby + ayd, + ay)2 
+ (ab; BO Wein ly to Gab AL andy + asb,)z~° (4.65) 
TO PICR eerie com eo co lesd, ce Dah ble Ne 
Chee + (ab, + a,jz + (aby + dgb, + a,b, + as)z~ 
4. (aib;, + a) Dye (tea 


—2 


Beginning with the output equation, equation (4.17), and referring to Figure 26, 
We note that the output of the multiplier d is dz-*X. This results is a revised output 


equation of 
Y=¢,4,+oV,+4zX. (4.66) 
Substituting the state equation into the output equation, results in 


¥ = [(b) + 303d, + Bec, + (b? + 2b,b,)c, Jz, 
+ [ (b)b, + 2b,bj)cy + (bp dy + BS)cy Jz V, 
+ oC). 
+ [agb, + ay)e, + ayia le ise 
aie [ (ayby + dgby + a,b; + ay)c, + (ab, + a)ep |z 7X 
+ [(agb + 2agbiby + ab} + aby + ayby)ey + (agbi + CN eT ae) (oral Fae @ 
ae [ (a,b; + 2a,6,6, + arb + a)b,)c, + (a,b Jo te ayby)cy 27 X 
a [ (ay; + 2ayb,by)c, + (a,b; + agby)c + aes 


(4.67) 


Y = [ (6; + 3252, + b3)c, + (by + 2b,b)c) Jz“ Y 
+ [ (bpd, + 2b,b3)c, + (bi by + Bey Jz Y 
+ ajc 
+ [layby + a )ey + acy ]27N 
+ [ (ab) + agb, + a,b, + ay)c, + (ap, + ay)ep ]J2-7-X 
as [ (ayb; + 2ajb, b> + a,b; + a,b, + a,b, )c, + (ab; + ab, + a,b, + een zea 


(4.68) 


+ | (a,b) + 2a,,b + agbj + andy), + (a,b] + ay by + ay, ey J27X 
+ [ (ab; + 2ayb,b,)c) + (a,b + Ayby)Cy + diate 


To demonstrate equality with the Y of the augmented IIR filter as expressed in equation 


ff oiore — 3, letc,=i,¢c,=0,andd=0. Then the output equation is: 


Y =(b, + 3b¢d, + be)z” YY + (bed, + 2b,bs)2 YY + aX + (qb, + a,)z X 
+ [ ag(b} +b) a0) + elem 
+ [ ag(b? + 2b,b.) + ay(b + b) + agd, J2z7X 
+ [a,(b} + 2b,by) + a(by + by) J27X + ay(by + 2b,b,)277°X. 


(4.69) 


By substituting: 


nN 
nN 


by ae Seah, oe = by, 
bib, + 2b,b3 = by 


Co Oi 
Agb, + a =a, 
ag(by + by) + ab) + a, = 43 Be) 


ag(b} + 2b,b,) + ay(by + by) + gb, = aya 
a,(b) + 2b,by) + a,(b? + by) = ay; 
a4(bj + 2b,b>) = ay¢ 


we can write this equation as: 


Y=5,,z Yt die 1 a, )X + aygz X 
ne ayy2¥ = ere ON (4.71) 
Bec te ayez X. 


This is the same as the equation (3.56) for p=3 on page 20. This realization is shown 
in Figure 26 on page 57. 

From equation (4.33) we can see that the feedback portion of the realization 
must have a delay of L stages in any loop. The general realization of that portion 1s 


shown in Figure 27 on page 38. 


FE. PIPELINED IMPLEMENTATION 

All of the IIR filters in section D can be developed into high-speed IIR filters by 
applying pipeline techniques. We use Ka-stage pipeline add units and Kn -stage pipeline 
multipliers to increase the speed of the filters. The Km stages are involved in the coefh- 
cient multiples and the Ka stages are involved in the adds. The net result is a (Km + 
Ka)-stage pipeline multiply-add unit. The equations for the filter outputs for different 
LS dre {1vem elon. 

lL L=1 


The equation is 


; oa i 12 | (4.72) 


+ a),X + ee Aone 


This reahzetic:ris shown in Figure 28 on page 59. 
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Figure 27. Heart of the Pipeline Recursive Filter Realization 


The equation 1s 
{—(2Ka+Km) > _-(2Ka+Km+1) 
zg ehatekm) = go kat2km) by\2 Ua by22 Me 
ia —ly 2An, —3 
+ Qi + Qppekek aes ee 2 ee 


This realization 1s shown in Figure 29 on page 60. 
3. L = 3 


The equation 1s 


= 0 F , . . , —, a = 
57 (l0Ka+2Km) y, _ ,—(10Ka+2Km) ca teemaee oa 2y 


=< yz X + Para ¢ 


This realization 1s shown in Figure 30 on page 61. 
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(4.73) 


(4.74) 





With Ka, Km 


] 


Second Order IIR Filter: L 


Figure 28. 





L = 2 with ka, Km 


Figure 29. The Augmented IIR Filter: 
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L = 3 with Ka, Km 


30. The Augmented IIR Filter: 


Figure 


61 


4 L=4 


The equation is 


bee nee bere ey, 


+ ana Ss Qa; ga te ayoz X 


This realization 1s shown in Figure 31 on page 63. 

From analysis of the feed back loop portion of the state space realization, as 
shown in Figure 27 and equation (4.33), we have seen that the delay around any loop 
must be L. A general realization for the feedback portion of the pipeline case 1s shown 
in Figure 32 on page 64. 


Therefore we can see that 
L=2Ka+t+ Km. (4.76) 


The following examples show the analvsis of the cases for L = 1, 2, 3 and 4. 
5. Example for L = 1: 
a ha=0,Am= 1 


L=2Ka+kKkm 
2(0) + 1 
=I 


I 


The only integer, non-negative values for Ka and Km which solve equation (4.76) are 
Ka=0,Am=1. This means we cannot use pipeline modules with longer delavs than 
Ka=0.Am=1. 

This realizes the original IIR filter for L = 1, but with an additional 2 de- 


lays and the equation 1s 
2° Ye= 27 [B20 YP Din? 0 eee nee cee (4.77) 


This realization 1s shown in Figure 33 on page 65. 
6. Example for L = 2: 
The onlv integer, non-negative values for Ka and Km which solve equation 
(4.76) are Ka = 0, Aim = 2 or Ka 
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Figure 


31. The Augmented IIR Filter; L = 4 with Ka, Km 


63 





Figure 32. Heart of More General Pipeline Filter Realization 


For these values, we use pipeline modules with delavs Ka=0,Am=2. 
This realizes the augmented IIR filter for L = 2, but with an additional 4 


delavs. The equation 1s 
ZY = z~ [bz °¥ + bam y+ Q,1X + day eek mee ee ances |. (4.78) 


This realization is shown in Figure 34 on page 66. 
b Aha=1,Am= 0 
L=2Ka+ Km 
= 2(1) +0 
=2 


For these values, we us pipeline modules with delays Ka=1, Km=0. 





l 


0, Km= 


IIR Filter: L = 1 with Ka= 


33. Second Order 


Figure 
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IR Filter: L = 2 with Ka=0, Km=2 


Figure 34. The Augmented I 


This realizes the augmented IIR filter for L = 2, but with an additional 9 


delavs. The equation is 
eee onze) ay rane + az. A ajgz (XX. (4.79) 


This realization is snown in Figure 35 on page 68. 
7. Example for L = 3: 
The onlv integer, non-negative values for Ka and Km which solve equation 
eo ate .a— 0, Ay =3 or Ka=1, km= 1. 
a ha = 0, Km = 3 


L=2Ka+Km 
= 2(0) +3 
=3 


For these values, we use pipeline modules with delavs Ka=0. Am=3. 
This realizes the augmented IIR filter for L = 3, but with an additional 6 


delavs. The equation is 


ae 
2 Sz 


(4.80) 


—3 r —4 Ld —| [2 
a6 6,42 i ta oes y+ Gy, x tn A122 x 
—2y- -3y, —4,- 
Se Qy3- x a Qyg< yy T Ay52 x 
This realization is shown in Figure 36 on page 69. 


b Kha= 1,Am= 1 


L=2Kat+ Am 
= 21) 4a 


sf 
2 


For these values, we use pipeline modules with delays Ka=1, Am=1. 


This realizes the augmented IIR filter for L = 3, but with an additional 12 


delays. The equation 1s 


—3 —4 J? a —| r 
ra =2 2 _ 2 , . 
aa Q)32 2h oe Qy4e ae oe 


This realization is shown in Figure 37 on page 70. 
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0 


2 with Ka= 1, Km= 


IIR Filter: L = 


The Augmented 


Figure 35. 





3 with Ka=1, Km=1- 


JIR Filter: L = 


The Augmented 


Figure 37. 


8. Example for L = 4: 
The only integer, non-negative values for Ka and Km which solve equation 
(4.76) are Ka =0, Am =4 and Ka=1, Km =2 or Ka =2, Km =0. 
a. Ka = 0,km = 4 


L=2Ka+ Km 
= 2(0)+ 4 
=4 


For these values, we use pipeline modules with delays Ka=0, Km=4. 


4, but with an additional 8 


This realizes the augmented IIR filter for L 
delays. The equation is 


-_ 


< 


—4 Ad —5 iz —]) ¥ 
a eee —§ by 42 )) a by22 Y+ ay 4X + A) 22 Xx 
i cae —2 -—3 yu —4 —5 T° 
=. Qy 32 X ar Qy4Z A 55 Qy5Z X+ Q162 x 
This realization 1s shown in Figure 38 on page 72. 
b. Kha = 1, km = 2 


L=2Ka+ Km 
= 2(1) +2 
=i 


For these values, we use pipeline modules with delavs Ka=1, Km=2. 
This realizes the augmented IIR filter for L = 4 but with an additional 15 


delays. The equation 1s 


—4 5, —5 —ly 
tae ris 63,2 } + 0452 Y+ ayX t+ A122 x 
° + ee ee aes eee ae: |X || 


This realization is shown in Figure 39 on page 73. 
c ha = 2, Km = 0 


L=2Ka+ Km 
= 2(2) +0 
= 4 


For these values, we use pipeline modules with delays Ka= 2, Km=0. 
This realizes the augmented IIR filter for L = 4, but with an additional 22 


delays. The equation 1s 


7) 





L = 





4 with Ka= 1, Km=2 


—— 


R Filter: 


The Augmented II 


Figure 39. 


—22 ? 
} 


Cyl biz POs) aero 
"| Mi 12 lel 12 (4.84) 


—?2 ry —3 J) —4 } A —5 oa 5 
7 Gig ox + Qy,4Z x + diez x + Ay 62 oe 


This realization is shown in Figure 40 on page 75. 


Table 5. STATE SPACE REALIZATION SUMMARY 


| iL | Ka-add units || Km-multipliers | D(output delayed) 
ME Gia | oe 31 


th2 | Gad |} Go 
f. | G2} bo 


ore 
Cat 


or) 





F. SUMMARY 

This chapter develops an alternate method for high speed implementation of state 
space digital filters. The proposed method implements A?’ instead of A as seen in 
equation (4.33). The proposed method allows arbitrary A in the implementation, and it 
can take full advantage of the optimum state space filters along with the benefit of A+ 


realization. The delay and complexity of each realization are summarized in Table 5. 
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VY. CONCEUSIONS 


A. SUMMARY 

This thesis describes pipeline realizations of the IIR notch filter. The transfer 
function and the coefficients that are appropiate for IIR notch filter have been pre- 
sented. Realization of a high speed IIR notch filter using a pipeline technique by Loomis 
and Sinha has been presented. An alternate technique in using state space realization 


by R. Gnanasekaran has also been presented. 


B. COMPARISON 

Table 4 in Chapter III summarizes the Loomis, Sinha pipeline recursive design and 
Table 5 in Chapter TV svuinmarizes Gnanasekaran’s state space realization. It was shown 
in Chapter III section I that we have only one choice of Ka and Km for each value of 
p and in Chapter IV section E that we have | choice of Ka and Km for L=1, 2 choices 
of Ka and Km for L=2 and for L=3, and 3 choices of Ka and Km for L=4. 

In the Loomius/Sinha pipeline recursive design we note that for p= 2, Ka=1, Km=2, 
output delav= 5, the filter is unstable. However. in Gnanaskaran’s state space realization 
for L=4, Ka=1. Km= 2, output delay= 15, the filter 1s stable. Thus we see that in cases 
where Loomis, Sinha realization are unstable, a state space realization always exists and 
is stable. However. it will require more hardware and have a greater overall delay. All 
of these results are summarized in Table 6. In that table, the column labelled D refers 
to the overall delay and the columns labelled #Mlt and #Add refer to the number of 
multipliers and number of adders as shown in the appropriate figure. 

Note that for a given value of Ka and km, the number of adders and multipliers 
required for the state space realization is larger than that required for the scalar. This 


disadvantage 1s offset bv the guaranteed stability of the state space realization. 
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Table 6. SUMMARY OF REALIZATIONS 


[ka [km [Lt [D | Fie [evn] ena] [D [Fie [em] ends 
Bo i 2 is ion |. TCT 
eet? |4 ie fs ts | [|| t. 
nfo f2 fs _f3 fis [3 fo [3 pi [3s 3 
ems (3 |¢ [sia lis | | (|. 1. 
nfs pe fs fi fas [i [4 [oo [es [e_ 
LC 
> fo [a (2 [ss po fr ||. |. [1 
Pes |. |. |. | is |e |» (: [s_ 
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APPENDIX A. PROGRAM AND DATA OF FREQUENCY RESPONSE 


A. PROGRAM FOR FREQUENCY RESPONSE CALCULATION 
DIGFREQ. FOR VERSION: 2/03/88 


PURPOSE: THIS PROGRAM COMPUTES THE FREQUENCY RESPONSE OF 
DISCRETE SYSTEMS. THE PROGRAM CONSISTS OF A MAIN 
PROGRAM THAT CONTROLS THE INPUT/OUTPUT AND THE 
SUBROUTINES dfresp AND coeff. SUBROUTINE dfresp 
COMPUTES THE FREQUENCY RESPONSE OF EACH SYSTEM. 
SUBROUTINE coeff ALLOWS THE USER THE OPTION OF 
GENERATING THE FILTER COEFFICIENTS OF THE SYSTEMS 
TO BE ANALYZED BY WRITING THE APPROPRIATE EQUATIONS. 

IF THE USER ELECTS TO GENERATE THE COEFFICIENTS BY 
USING THE SUBROUTINE coeff, THE EQUATIONS MUST BE 
WRITTEN INTO THE SUBROUTINE USING STANDARD FORTRAN 77 
STATEMENTS. THE COEFFICIENTS MUST BE STORED IN THE 
ARRAYS b() AND c() WHICH CORRESPOND RESPECTIVELY TO THE 
NUMERATOR AND DENOMINATOR TERMS OF THE SYSTEM EQUATION. 
THE USER CAN SELECT ONE OF TWO OPERATING MODES: BATCH 
OR TEST. IN BATCH MODE THE AMOUNT OF INTERFACE WITH 
THE USER IS MINIMIZED AND IT IS ASSUMED THAT THE INPUT 
DATA HAS BEEN STORED IN THE DEFAULT FILE 'DIGFREQ. IN’. 
IN TEST MODE THE USER IS PROMPTED FOR THE NAME OF THE 
INPUT FILE OR HAS THE OPTION TO PERFORM A TRIAL RUN BY 
USING THE INPUT DATA STORED IN THE FILE 'DIGFREQ. TST’. 
IT IS RECOMMENDED THAT FIRST-TIME USERS SELECT THE TEST 
MODE AND MAKE A TRIAL RUN WITH THE PRESTORED INPUT DATA. 
THE TEST MODE ECHOES THE INPUT DATA ONTO THE MONITOR TO 
ALLOW VERIFICATION OF ITS ACCURACY. THIS PROGRAM WILL 
COMPUTE THE FREQUENCY RESPONSE OF UP TO THREE SYSTEMS. 
FOR EACH SYSTEM, THE USER HAS THE OPTION OF HAVING THE 
OUTPUT EXPRESSED IN DECIBELS (db). THE OUTPUT OF THIS 
PROGRAM IS STORED IN TABULAR FORM IN THE FILE 
"DIGFREQ. OUT' AND IN A FORM SUITABLE FOR PLOTTING 

IN THE FILE 'DIGFREQ. DAT’. 


WETTED TERETE TEE TES TN PUT  vededededodevededevededodivevedevere devon veve ere rere e ee 


THIS PROGRAM ASSLMES THAT EACH DISCRETE SYSTEM IS MODELED BY THE 
EQUATION: H(z) = num/den WHERE: 


num = b(0)*z**1, + pC 1)*2**( ba] et B=) 2 ecto 


den = c(0)*z**N + cC1)*z**(N-1) +... + cCN-1)*z + c(N) 

L = A NON-NEGATIVE INTEGER, THE DEGREE OF THE NUMERATOR 
POLYNOMIAL. 

A NON-NEGATIVE INTEGER, THE DEGREE OF THE DENOMINATOR 


Ca Gurca Ga ©) GauGa CG) Ca. Coie) Gs Cy C2 Ca C2 Co GO) Ca Co ae? Ci) Ca Gr Ca GC C7 ee) Os CG) Ga) Gt 
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POLYNOMIAL. 
b(0)...b(L) = REAL COEFFICIENTS OF THE NUMERATOR TERMS. 
c(0)...c(N) = REAL COEFFICIENTS OF THE DENOMINATOR TERMS. 


THE INPUT PARAMETERS SHOULD BE STORED IN A FILE NAMED 
'DIGFREQ. IN’. ALL OF THE READ STATEMENTS USED BY THIS PROGRAM 
REQUIRE FORMATTED INPUT. PARTICULAR ATTENTION SHOULD BE PAID 
TO THE FORMATS, ESPECIALLY THE USE OF THE DECIMAL POINT TO 
DENOTE 'REAL' NUMBERS. THE INPUT PARAMETERS REQUIRED BY THE 
PROGRAM ARE LISTED BELOW. 


NAME yes RANGE (ARRAYS) RESTRICTIONS 
numsys INTEGER 1 <= numsys <= 3 
L INTEGER 0 <= L <= 128 
N INTEGER 0 <= N <= 128 
dsorce CHARACTER = Oke. 
yscal CHARACTER Poll POR” Hoge 
theta0 REAL 

dithta REAL 

numpts INTEGER 1 <= numpts <= 101 
b() REAL OF le 2e » L 0 <= L <= 128 
et) REAL OC; 1), gan aN 0 <= N <= 128 
WHERE: 


numsys = THE NUMBER OF DISTINCT SYSTEMS H(z) TO BE ANALYZED. 
tiie INTEGER VALUE MUSE OGGUR AleaHE TOP OF THE INPUT 
Pere) DELINEATES tHE NUMBER OP SYSTEMS TO BE READ BY 
THE PROGRAM AND ANALYZED. FOR EACH SYSTEM (1, ..., numsys) 
iPr onaibleERS BELOW MUST APPEAR IN THE INPUT FILE. 


L = AN INTEGER VALUE SPECIFYING THE DEGREE OF THE NUMERATOR 
POLYNOMIAL. 
N = AN INTEGER VALUE SPECIFYING THE DEGREE OF THE DENOMINATOR 


POLYNOMIAL. 
Gsorce = THE CHARACTER STRING 'F' OR 'S'’ DENOTING WHETHER THE 
SYSTEM COEFFICIENTS ARE TO BE READ FROM THE INPUT FILE (F) 
OR GENERATED (S) THROUGH USE OF THE SUBROUTINE coeff. 
yscal = A CHARACTER STRING SPECIFYING THE DESIRED MAGNITUDE OPTION: 
"STD' WILL PRODUCE STANDARD MAGNITUDE OUTPUT; 
'LOG' WILL PRODUCE MAGNITUDE EXPRESSED IN DECIBELS (db). 
thetaO = THE STARTING VALUE OF THETA (RADS) AS IN Z=EXP(J*THETA). 


dltnta 


Il 


THE INCREMENT OF THETA (RADIANS). 


numpts = THE NUMBER OF FREQUENCY POINTS FOR WHICH THE OUTPUT IS 
TO BE COMPUTED. 


BO ee Ieee RATOR@COLEFICTENDS IN ORDER b(@), b(1), ..., BCE): 


IF dsorce = 'F' IS SELECTED THEN THE USER MUST SUPPLY THE 
L+1 NUMERATOR COEFFICIENTS IN THE FILE. IF dsorce = 'S' 


co 


Corr Gel Gela Gel Ge GueGaN Gl Gung Gl ck Ga 1G) Ga Gl ca to -GacG Ca Ga Gl CavC) G) °C) Cac) CG) Co Ca C2) G) Ca-Ca Cl CG) Ca: Ga GaiGarga cd GajGasG. Ca GicG) Cisco Garc arc) 


THEN THE USER HAS ELECTED TO GENERATE THE NUME Rater 
COEFFICIENTS BY WRITING THE APPROPRIATE FORTRAN STATEMENTS 
IN THE SPACE PROVIDED IN SUBROUTINE coert: shi eiii se tHOD 
OF DATA GENERATION IS ELECTED THE PROGRAM MUST BE RECOMPILED 
BEFORE EXECUTION. 


c() = THE DENOMINATOR COEFFICIENTS IN ORDER ec (0) cs) sere One 
IF dsorce = 'F' IS SELECTED THEN THE USER MUST SUPPLY THE 
N+1 DENOMINATOR COEFFICIENTS IN THE FILE. IF dsorce = 'S' 
THEN THE USER HAS ELECTED TO GENERATE THE DENOMINATOR 
COEFFICIENTS BY WRITING THE APPROPRIATE FORTRAN STATEMENTS 
IN THE SPACE PROVIDED IN SUBROUTINE coeff. IF THIS METHOD 
OF DATA GENERATION IS ELECTED THE PROGRAM MUST BE RECOMPILED 
BEFORE EXECUTION. 


NOTE: THE INPUT FORMAT STATEMENTS OCCUR IN THE MAIN PROGRAM 
FOLLOWING THE CAPTION: *etdtretordee = =TNPUT FORMAT = *icteotvereaeet, 
THE FORM OF THE INPUT DATA FILE IS: 








LINE # ENTRIES FORMAT 
1 numsys vee | 
2 L,N,dsorce,yscal 13,01 Peer al jt3 las 
3 dithta,theta0,numpts Z2eu0nO 513 

Nod ie bCk).k=Onige ee Lb 6f 1102 0 

NOE 2 CUK)ey k= 05 eee Snel e 

NOTE 3 

WHERE: NN + (L/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER). 


ND + (N/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER). 
NOTES 1. THE NEXT NN LINES ARE ONLY REQUIRED IF dsorce = 'F'. IF 
dsorce = 'S' THEN THE USER HAS ELECTED TO GENERATE THE 
Lt+1 NUMERATOR COEFFICIENTS IN THE SUBROUTINE coefé. 
THE USER MUST PROVIDE THE APPROPRIATE FORTRAN STATEMENTS 
IN SUBROUTINE coeff TO GENERATE THE VALUES FOR b(). 


2. THE NEXT ND LINES ARE ONLY REQUIRED IF dsorce = 'F'. IF 
dsorce = 'S' THEN THE USER HAS ELECTED TO GENERATE THE 
N+l DENOMINATOR COEFFICIENTS IN THE SUBROUTINE coeff. 
THE USER MUST PROVIDE THE APPROPRIATE FORTRAN STATEMENTS 
IN SUBROUTINE coeff TO GENERATE THE VALUES FOR c(). 


3. FOR numsys > 1 THESPORMAT ORSUIND Se aloes beach). 


4. THE FORMAT £10.0 USED FOR INPUT DATA PERMITS THE DECIMAL 
POINT TO BE PLACED ANYWHERE IN THE FIELD OF 10 COLUMNS 
AND ALSO ALLOWS THE EXPONENTIAL FORMAT TO BE USED (CEG. 
3146.2 = 3. 1462E+03). 


VEDETE TE TODD ETTORE TETRA OUTPUT — vedededededevedodevedevededede done tecededede de tenet enone ae 


THE OUTPUT DATA CREATED BY THE PROGRAM IS STORED IN TABULAR FORM 
IN THE FILE 'DIGFREQ. OUT’. ADDITIONALLY, THE OUTPUT DATA IS 
WRITTEN INTO THE FILE 'DIGFREQ. DAT' TO FACILITATE PLOTTING BY 
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A SEPARATE, USER SUPPLIED PROGRAM. THE FORMAT OF THE DATA IN 
'DIGFREQ.DAT' IS: e12.6, 2x, e12.6. THE FIRST ENTRY CORRESPONDS 
TO THE ORDINATE VALUE (THETA) AND THE SECOND ENTRY THE ABSCISSA 
VALUE (MAGNITUDE OR PHASE). ADDITIONAL HEADER INFORMATION IS 
WRITTEN INTO THE DATA FILE TO ALLOW FOR CONTROL AND LABELING OF 
BAGh PLOT. 


THE INPUT PARAMETERS FOR THE SYSTEM DESCRIBED BELOW ARE STORED IN 
THE SAMPLE INPUT FILE 'DIGFREQ.TST' AND CAN BE USED FOR A TRIAL 
RUN IN THE TEST MODE. 


SYSTEM: H(z)=z/(z-0.5) 

GOAL: TO OBTAIN THE FREQUENCY RESPONSE FOR THIS SYSTEM FROM 
THETA = 0.0 TO THETA = 3.14159 (PI RADIANS) IN STEPS 
OF dithta = PI/10.0 


Bele tHE SsYolEt! DESCRIBED ABOVE THE INPUT FILE IS: 


1 

001 001 F STD 
ge oi1g159° "0. 0 Ort 

0 0.0 

da 8 =O.) 


THE RESULTING OUTPUT DATA FILE: 'DIGFREQ. OUT' IS: 
INPUT DATA FOR SYSTEM # 1 


INPUT DATA SOURCEFILE: DIGFREQ. TST 


DEGREE OF NUMERATOR = il 

DEGREE OF DENOMINATOR = 1 

dsorce = F 

NUMBER OF FREQUENCY POINTS = li MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA =~ . O0QOQOOQ00Q0E+09 

INCREMENT OF THETA = .314159E+00 


THE NUMERATOR COEFFICIENTS b(0),b(1)...b(L) ARE 


mOOOr-O1 . OOOOE+00 


THE DENOMINATOR COEFFICIENTS c(0),c(1)...c(N) ARE 


. LOOOE+01 =, DOOCETOO 
OUTPUT DATA FOR SYSTEM # 1 


THETA MAGNITUDE PHASE 
(RADIANS ) (DEGREES ) 


8 | 


. OOOO00E+00 . 200000E+01 . OOOOO0E+00 
- 314159E+00 8269 7E0 | . 164149E+02 
. 6283 18E 00 - 150588E+01 - 260267 /EOz 
.942477E+00 PI 2Z7660E 70" » 29607 Ea02 

. 125664E+01 . LO3088E+01 . 293546E+02 
. 15 7060E401 - 894428E+00 , 26565 1b+02 
. 188495E+01 . 800894E+00 - 223602E 102 
2IPDLIE Eon . 737654E+00 / 17 3606E+02 
25 le2/ Ee . 696900E+00 . 118186E+02 
2027 4aSb rol . 674038E+00 09772 3b-01 
- 314159E+01 . 666667E+00 -484184E-04 


eee END OF RUN, SYSTEM #1 ------------- 


FOR ILLUSTRATIVE PURPOSES THE COEFFICIENTS b() AND c() COULD 

HAVE BEEN GENERATED BY SPECIFYING dsorce = 'S' AND WRITING THE 
APPROPRIATE FORTRAN STATEMENTS INTO SUBROUTINE coeff. THE 
STATEMENTS THAT COULD BE USED TO ACCOMPLISH THIS ARE WRITTEN INTO 
THE SUBROUTINE BUT ARE 'COMMENTED OUT’. 


C2 GaGa Oo Ca GC) Ca) C26) GC): C) Cle a Gs) GGG) 


vervededederedesodededededesee te seteseiesevee ae = MATN PROGRAM  witdededededsosiededtatdsdedededevededevedesedede dea 


character infile*12, mode*1, ylabi*13, dsorce*1, yscal*3 
real mh(101), ph( 101), thetav(101), c(0:128), b(0: 128) 


C PROMPT USER FOR MODE: BATCH OR TEST. 


write(*s 1 Pils) 
read(*,1117) mode 
if((mode. eq. 'Y'). or. (mode. eq. 'y')) then 
mode =  y" 
write(*,1118) 
read(*,1119) infile 
else 
infiie = 'DIGFREQ. IN' 
endif 


C UNIT=1 DEFINED AS INPUT FILE. UNITS=2735 DEFINED As OUT Or rire 


open(unit=1,file=infile,status='old' ,iostat=ierr, err=999) 
open(unit=2,file='DIGFREQ. OUT' ) 
open(unit=3,file='DIGFREQ. DAT' ) 


C READ INPUT PARAMETERS AND CONDUCT ERROR CHECKS. 


read(1,1000) numsys 
numplts = numsys*2 
write(3,2000) numplts 
if((numsys. lt. 1). or. (numsys. gt. 3)) then 
write(*,1122) numsys 
stop ‘Error, numsys must be in the range: 1 <= numsys <= 3. 
endif 





do 10 nsys=1, numsys 
data mh/101*0.0/, ph/101*0.0/, thetav/101*0. 0/ 
ama, 1290.07, c/129*0. 0/ 


nread(1,1001) L, N, dsorce, yscal 
read(1,1002) dithta, theta0, numpts 


ieee 0), or. (L.et.128)) then 
write(*, 1124) nsys, L 
eGpmenLron, Limust be in the range: 0 <= L <= 128. ' 
eltecan CN. lt. 0). Orv CN, ct.1259) eehen 
write(*, 259) 1S ysiguN 
stop ‘Error, N must be in the range: 0 <= N <= 128. 
endif 


if((dsorce. eq. Po menGdsereescas £ ))) then 


fiserce = fF’ 
elseif((dsorce. eq. §'). or. (dsorce. eq. 's')) then 
isorce — S| 
else 
merece. *. 1018) crorce 
Seem tiesallowed values for dsorce are: ‘'S'' or ‘'F’'.' 
endif 


feeenumpts. it. 1).or. (numpts. gt. 101)) then 
write(*, 1127) nsys, numpts 
stop ‘Error, numpts must be in the range: 1 <= numpts <= 101.' 
endif 


if((yscal. eq. “STD! ). or. (yscal.eq. 'std')) then 
De 


yscal = ‘ST 
ylabl = ' MAGNITUDE ' 
elseif((ys cal. eq. LOG )wor.(yseal.,eq) loge ))ethen 

yscal = ‘LOG’ 

ylabl = "MAGNITUDE( db) ' 
else 

write(*,1128) yscal 

MeOpmetELOL, yscal must be the string: ‘L0G ' or STD 9. 
endif 


C FOR dsorce = ‘'F' READ THE COEFFICIENTS b() AND c() FROM THE INPUT 
C FILE. FOR dsorce = 'S' CALL coeff TO GENERATE THE COEFFICIENTS. 


if(dsorce.eq. 'F') then 
read(1,1003) (b(k),k=0,L) 
read(1,1003) (c(k),k=0,N) 
else 
call coeff(L,N,nsys,b,c) 
endif 


Coevelne INPUT DATA INTO THE OUTPUT FILE: DIGFREQ:. OUT. 
write(2,1008) nsys 
write(2,1010) infile 


werte( 2,1110) L 
write(2,1111) N 


§3 


Write(2, 1019 edserce 
write(2,1112) numpts, yscal 
write(2,1113) theta0d 
write(2,1114) dlithta 
write(2,1004) 
write(2,1005) (b(k),k=0,L) 
write(2,100G) 
write(2,1005) (c(k),k=0,N) 
write(2,1009) nsys 
write(2,1126) ylabl 
write(2,1007) 


C FOR TEST MODE ECHO ALL INPUTS ONTO MONITOR (UNIT = *). 


if(mode.eq. 'Y') then 
write(*,1120) nsys, infile 
write(*,1110) L 
Wirieel-s iolela)igy 
writec*,1019) dsorce 
write(*,1112) numpts, yscal 
write(*,1113) thetad 
write(*,1114) dlthta 
write(*,1004) 
write(*,1005) (b(k),k=0,L) 
write(*,1006) 
write(*,1005) (c(Ck),k=0,N) 
write(*,1123) nsys 
pause 'END OF RUN. STRIKE <CR> WHEN READY TO CONTINUE. ' 

endif 


C CALL dfresp TO CO“DUTE THE FREQUENCY RESPONSE. 
call dfresp(b,c,mh,ph,L,N,theta0,dlithta,thetav,numpts,yscal) 
© WRITE RESULTS INTO VOULPUD FILE: SOvGrREG ar: 


write(3,2001) numpts 
write(3,*) ‘MAGNITUDE RESPONSE’ 
write(3.:) me SHELAe CR sDeye 
write(3,2003) ylabl 
do 55 np=1, numpts 

write(3,2010) thetav(np), mh(np) 
55 continue 


write(3,2001) numpts 
write(3,*) ‘PHASE RESPONSE’ 
write(3,*) ‘THETA (RADS)' 
write(3,2003) ' PHASE (DEG) ' 
do 56 np=l1, numpts 

write(3,2010) thetav(np), ph(np) 
56 continue 


GC WRITE RESULTS INTO OUTPUT BILE DICGrREO Oui 


do 150 np=1, numpts 
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write(2,1013) thetav(np), mh(np), ph(np) 
SO continue 
write(2,1123) nsys 


LO continue 


write(*,1121) 
999 close(unit=1) 
close(unit=2) 
close(unit=3) 


gi(ierr, st.0) then 
write(*,1116) ierr 
encif 


Cisveisiede: seis INPUT FORMAT creel ees soy aby sriaeisriac bariovlariselee ts 


1000 format(il) 

Be Geet ormat( i3,tl1,i3,t21,a1,t31,a3) 
mooze rormat( 2£10.0,i13) 

1003 format(6f10. 0) 


(Okikitikitikieteititistieisisis eb tite hiss tl sei dle el eit flee el Slee eels 


1004 format(t4,' THE NUMERATOR COEFFICIENTS b(0),b(1)...b(L) ARE: 

emo tormat(6(2X,e11.4),//) 

1006 Oformat(//,t4, ‘THE DENOMINATOR COEFFICIENTS c(0),c(1)...c(N)' 
rc ARE /) 

1007 format(t6,' (RADIANS)',t38,'(DEGREES)',/) 

1008 format(t16,' INPUT DATA FOR SYSTEM # ',i1,//) 

1009 format(///,t16,' OUTPUT DATA FOR SYSTEM # ',i1,/) 

1010 format(t4,' INPUT DATA SOURCEFILE: ',al2) 

Hols format(t4,3(e12.6,4x)) 


iMetees format( ix, dsorce = ',al,2x, Error, illegal value for dsorce. ' 


1019 format(t4,'dsorce = ',al) 

1110 format(t4,' DEGREE OF NUMERATOR = ',i3) 

1111 format(t4,'DEGREE OF DENOMINATOR = ' ,i3) 

1112 Oformat(t4,'NU:°3ER OF FREQUENCY POINTS = ',i3,t39,'MAGNITUDE' , 
PormION = as) 

imtoo format(t4,’ STARTING VALUE OF THETA = ' ,e12. 6) 

1114 format(t4,' INCREMENT OF THETA = ',e12.6,/) 

1115 Oformat(1x,' vO YOU WISH TO RUN THIS PROGRAM IN TEST’, 
1' MCDE ? GN ye<CR> 2 2) 

1116 Oformat(///,1x,' ERROR OPENING INPUT FILE, PROGRAM TERMINATED. ' 
Weis ERROR CODE: ,i4,/////) 

ily SStormat( al) 

1118 Oformat(/////,1x,'TYPE THE NAME OF YOUR DATA FILE FOLLOWED’ , 
1' BY <CR>.',/,' IF YOU DESIRE TO MAKE A TEST RUN USING THE’, 
2' SAMPLE DATA ALREADY STORED’ ,/,' IN THE FILE: DIGFREQ.TST' , 


; ivERoePIGEREOG. 25m <CR>../. FILENAME: °, ,) 
1119 format(al2) 
1120 format(////,t4,' SYSTEM samc ile INPUT DATA SOURCEFILE: ', 


TUE Oformat(//, t4,'TABULAR OUTPUT DATA IS STORED IN FILE: DIGFREQ.OUT', 


lve + ' PLOTTING DATA lSeolORED IN FILE: DIGFREQ. DAT. ) 
22 ene) //. pee nes valivesotunumsys iss 9.11, . ») 
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pss 
1124 


iS 


ize 
a2 
LEZs 
2000 
2001 
2003 
FAVA, 


a TS Sa 1s Jd Re 


format(/,1x,13('-'),' END OF RUN, SYSTEM =) saemeoaelc (gees {/) 


Oformat(//////, Ge ‘The degree(L) ae the numerator for System , 


ler eile 1S 6 it = cts ay 
Oformat(//////, 2 ne “degree(N) of the denominator for system’ 
Le gh ae aS ee ON see 


') 
format()/ 7) sco. ‘THETA’ ,t21 al3, t40,'PHASE') 
format(/////7 ten lhe “value of numpts for tas | cl, See Stee aon) 
format(/////,t25 lhe value of yocalerce yao 
format(il) 
f£Ornmatt it) 
format(al3) 
format(el2.6,2x,e12. 6) 


end 


SUBROUTINE: dfresp 


PURPOSE: THIS SUBROUTINE COMPUTES THE FREQUENCY RESPONSE OF 


THE SYSTEM. ALL FREQUENCY CALCULATIONS ARE IN RADIANS, 
HOWEVER THE OUTPUT IS CONVERTED TO DEGREES. 

THE OUTPUT FORMAT FOR EACH FREQUENCY INCREMENT IS: 
MAGNITUDE(M) PHASECE) AS IN: @aeEXP(J*P). 


subroutine dfresp(b,c,mh,ph,L,N,theta0,dlthta,thetav,numpts,yscal) 
reel mh(numpts), ph(numpts), thetav(numpts), imz, rez 

mea le biGOt)..-CCOaN) 

character yscai*3 

complex z, den, num, h, ci 


G DEFINE CONSTANTS. 


jeu 


4. O*%atan(1.0) 
CimOmn. 0) 


raul 


GC ITERATE FROM theta0, IN INCREMENTS OF dIthta. 


num 


den = 


do 100 np=1, numpts 
ci*b(0) 
enlcore (0) 


thetav(np) = thetad + (np-1)*dlthta 
rez = cos(thetav(np)) 


imz 


Sin(thetav(np)) 


Z = emplxCrez, 1mz) 


C CALCULATE NUMERATOR FOR GIVEN VALUE OF THETA, IF L > 0. 


if(L. gt.0) then 
do. 50°k=1.21 


50 


num = z*num + cip(k) 


continue 


endif 
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C CALCULATE DENOMINATOR FOR GIVEN VALUE OF THETA, IF N > 0. 


Mmenece. 0) then 
do 70 k=1, N 
den = z*den + ci*c(k) 
70 continue 
endif 
h = num/den 


C CONVERT COMPLEX VALUE 'h' INTO MAGNITUDE(mh) AND PHASE(ph) TERMS. 
C IF yscal = 'LOG' THEN CONVERT MAGNITUDE TO DECIBELS (db). 
C DIVIDE BY ZERO AVOIDED BY 'if' STATEMENTS. 


mh(np) = cabs(h) 
if(yscal.eq. 'LOG') then 
if(mh(np). gt.0.00001) then 

mh(np) = 20. 0*logl0(mh(np)) 
else 
mh(np) = -100.0 
endif 
endif 


if(abs(real(h)).1t.1.0e-15) then 
if(abs(aimag(h)).le.1.0e-15) ph(np)=0. 0 
if(aimag(h). gt.1.0e-15) ph(np)=90. 0 
if(aimag(h).1t.-1.0e-15) ph(np)=-90. 0 
else 
ph(np) = (180. G/pi)*atan2( aimag(h),real(h)) 
endif 


HOO CONntinue 


Tecuyrn 
end 


C SUBROUTINE: coeff 


PURPOSE: THIS SUBROUTINE ALLOWS THE USER TO GENERATE THE 
NUMERATOR AND DENOMINATOR COEFFICIENTS THAT DESCRIBE 
EACH SYSTEM TO BE ANALYZED. IF dsorce = 'S' THEN 
THE MAIN PROGRAM WILL CALL THIS SUBROUTINE. 


C7? a ko Ca 


Stproutine coefr(L,N,nsys,b,c) 
meal (0: L), c(0:N) 


pi = 4. 0*atan(1l. 0) 
el =L 
en = N 


owen aulentauanteutsciauaseselsuiautanaen'sulawaulsu's Pe Sawawton PP fa ale a cd - = = - a ate -f et = 
Crstevededevetete dete sete terete le Tete cede TOPE LETTS TS TCE TET OTS TC TE TS ETE TE TCT TET TER TOTES TCE TET TET ETE ETE TCT. 


C DEVELOP THE EQUATIONS TO GENERATE VALUES FOR THE ARRAYS b() AND c() 
C IN THIS SPACE. THE STATEMENTS TYPED IN MUST FOLLOW STANDARD 
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fh 


C2 C2 Cece) Co CCC), CG cn) Cac) 
Loe) 


FORTRAN 77 RULES AND MAY USE FORTRAN 77 INTRINSIC FUNCTIONS SUCH AS: 
SINC), COSC), ABS(C)... AN EXAMPLE IS SHOWN BELOW. NOTE THAT THE 
VALUE nsys CAN BE USED TO DISTINGUISH BETWEEN SYSTEMS IF MORE THAN 
ONE SYSTEM (numsys > 1) IS TO BE ANALYZED. 


vedeve = EXAMPLE acacac 


if(nsys.eq.1) then 
do 2 i=0, L 
b(i) = cos(i*pi/(2.0*el)) 
continue 
do 3 i=0, N 
c(i) = cos(2. 0*1i*p1/@o,0*en) } 
continue 
endif 


Cae rere seve de sede se dese se Fede He TEFEN Te eT ETT TENTS TEN FETE HE TET To To TEE BIE TERE TCR TRE TK EEE EEE EE BEER REE 


|o%, 


return 
end 


MAGNITUDE, PHASE OF SECOND ORDER IR NOTCH FILTER: P=0 


INPUT DATA FOR SYSTEM # 1 


INPUT DATA SOURCEFILE: digfreq.tst 

DEGREE OF NUMERATOR = 2 

DEGREE OF DENOMINATOR = 2 

dsorce =F 

NUMBER OF FREQUENCY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = .OQOQO0O00E+00 

INCREMENT OF THETA = _ .280600E-01 


THE NUMERATOR COEFFICIENTS b(0),b(1)...b(L) ARE: 


- 9900E+00 =7o209E700 27 OOE+00 


THE DENOMINATOR COEFFICIENTS c(0),c(1)...c(N) ARE: 


-LOOOE-OT So EO U »2801E+00 


8§ 


THETA 
(RADIANS) 


. OOOOQ0E+00 
e250000E-01 
ete OOR-O1 
. 841800E-01 
et 2 240E4-00 
male OS 00E4-00 
mes SO0E+00 
. 196420E+00 
. 224480E+00 
m2 2 5+0E+00 
-. 280600E+00 
. 308660E+00 
- 33672Z0E+00 
. 364780E+00 
23926405100 
- 420900E+00 
- 448960E+00 
7477 OZ0E+00 
=o UD 060E700 
.533140E+00 
pool Z00E+T00 
#2892Z60E700 
mol 73 Z0E+00 
- 6453580E+00 
. 673440E+00 
~ 7 01500E+00 
m7 22 300E +00 
Hi ST O2UETOO 
mos oc UrsO00 
Hoes? 4 OE+O00 
- 841800E+09 
. 869860E+00 
262722051700 
25 OUEs-00 
- 954040E+00 
- 982100E+00 
momo ToOE+O1 
oie eZ2E+0 1 
~eeooZ8E+01 
. 109434E+01 
.112240E£+01 
SiS O02oE+01 
see GoZh+01 
~iavosorns0 1 
. 123465E+01 
pe 27 OE 1 
mee 7 OC 1 


MAGNITUDE 


2999939 TO o 
£9999 39 Eo 
59999 39 Ene 
S999 75 7Rado 
. S992 SOE Oo 
-99 97 SSErOU 
s99DO STEAD 
<999936E 400 
29999550500 
. 999934E+00 
, 99993 3E70C 
; OOS OS TEs 
, 999930Es 09 
999978 T00 
5 ee She elo ONG, 
 IIDI Ze raeOe 
woe 22 Une) 
P99 Eee 
. 9999 lareoo 
- 9999 10E 08 
; 9999062706 
: GO9DOZE Aas 
- 999 BV CE UG 
-999S It eGo 
. 999884E+00 
99967 ORO 
- 99S8oSr00 
, 99S S5E4o0 
. 999847E+00 
, 999CS5E G00 
- 9973208700 
. 999S0SEa00 
9997 3b 00 
. 999/59 Ea 
«99373157500 
. 999697 ETO0 
«99905 ou Ro0 
. 999603ET00 
~999537E R00 
- 999453E+00 
. 999343E700 
H9SI19 SEU 
29 9b 65 Ea00 
. 92668 1E+00 
e727 621 3E700 
. 997440E+00 
. 996024E+00 


OUTPUT DATA FOR SYSTEM # 1 


PHASE 
(DEGREES ) 


- QOOO00E+00 


89 


219394 /E=0) 
w2OC2G 7c 01 
. 983436E-01 
-7/97845E-01 
»977/54h 201 
117779 E209 
2 1SC0S0E IG 
S1S8>6TETeO 
. 179474E+00 
e20075GE 400 
w2225 1604.00 
. 244780E+00 
9207020800 
e271 1IZE CG 
Mol SoZ Eee 
. 340366E+00 
. 366300E+00 
7393251 ERoD 
- 421327E+00 
. 450663E+00 
.481405E+00 
P2137 23a Oo 
- 547804E+00 
-o O56 72200 
7022182600 
sicher hls Secs 818, 
. /06764E+00 
1537 SEO 
- 804587E+00 
(O59 743Er 00 
- 919964E+00 
77501 00E.00 
LUD o 2 LEOt 
11 406HE 01 
PZ Sie eee 
pis 3oSCh od 
S145 309E io 
719909 SE ao | 
el7S176E01 
. 194344E+01 
7216202 t0m 
. 246560E+01 
mae Gye) dace 
“33 Z65eEtO0l 
. 400969E+01 
POU2Soon01 


to NSS 7 Ou 7 9IZ99SE TOO =0 7424 Eo 
»1346SSE+01 . 984530E+00 “21 00259h402 
- 137494E+01 - 942857E+00 = 194059E02 
; 140300E ea youlIUSE-Ol =O7 JOO JE nU2 
. 143106E+01 797679 20E700 . 200014E+02 
eke oo alee Oa . 984111E+00 plOZ595E+O2 
. 148718E+01 pO Z 607 EROO NOS 7 OLE SOd 
Peles ALS 27132210) 11 p99597 TE Ree - 1611S reeal 
, 1943308 40m 79974 SE Ree .413001E+01 
Res ile: o 210i - IIS 1IISEROG . 344006E+01 
. 1599428401 - 99367 250° 2945468401 
. 162748E+01 99897 2 Oe 125 / S10RaCal 
, 165554E 701 . 999190E+00 m2 25226) 50m 
; 168360R701 . 999341E+00 . 204852E+01 
Bey lilicielae a0) - 99945ZE*00 7 165629E208 
P97 2 tees Ol 999537 Ese l69SZ21EeO 
ly O77 ob eda . 999603200 p155806E Raa 
17 958ere 00 . 999655ET00 - 143975E-01 
, 18 2390R 0! - I996ISE ROG ; 133653ER01 
,16519Cr rod - 999 /SZE Ue .124550E+01 
. 188002E+01 , 999761 Ew . 116455E+01 
LPOG 0SEs01 » 9997 65E08 2 109 1 96E a Gal 
1936142 0n - 999805E500 . 102647E+01 
1964208401: -999822E noo . 966960E+00 
S98 22e8s0i . 99983 7E E00 : 91257 9R 60 
5 0 0 lan O)a . J99S650E Oe - 6626328700 
Pou eooGE OM . 999 861E Oo »816531E+00 
. 207644E+01 . 99967 IEoe p77 37 92E Ow 
. 210450E+01 - 999SS0E Ue . 734015E+00 
o21g2ocErOt - JIDSSOEEOU . 696843E+00 
peEOOUO2E EO - 999894E+00 ©001967E+00 
2165608701 : 999 EE eo mo2919SELOo 
© Ea Moeeie 10) - 9999 06ET0G . 598246E+00 
. 224480E+01 - 9999 ViEEee ele OeoU0 
ee 2OCn Ul 9999 >E Ree .541144E+00 
Peo O09ZErOM 9999 19 EOO . 514683E+00 
»2320708 70H : 9999 23r new .489434E+00 
» 235704E+01 . 99992Z6E TOO .465285E+00 
ee SSO IT a0) ~999929R Oo .442146E+00 
. 241316E+01 29999325 Oo P4199 Teh ao 
244 122E50u : I99DSSGETOG 7 O9652 12 Ow 
=240928E+01 1 9907o7 LOO 1 3//7689E 200 
B2o97 54810) 999929 EBOe 2357941800 
122550801 - 999941E+00 yo [OO SI ERO” 
he /aore wer (9) 050), . 999942 E Oe S199 ISE OO 
pec en OM . 999944E+00 - 301724E+00 
. 260958E+01 - 999945800 . 284024E+00 
S205 FOCHAO1 - 999947 E00 3206/77 1800 
BZ O65) OF FO: - 999945500 » 249979 n aw 
eos 7 OL FOI - 999949E+00 . 233459E+00 
ee? 2 oe el - 999950R R00 .21/33ZER 09 
O27 4756E.0n - 9999510500 »201519E aCe 
2/7 7 OGETO - 9999525509 7, 1859878700 


END OF RUN, SYSTEM # 1 


90 


C. MAGNITUDE, PHASE OF THE AUGMENTED IIR NOTCH FILTER: P=1 


INPUT DATA FOR SYSTEM # 1 


INPUT DATA SOURCEFILE: digfreq.tst 
DEGREE OF NUMERATOR = = 3 
DEGREE OF DENOMINATOR = 3 
dsorce = F 

NUMBER OF FREQUENCY POINTS = 
STARTING VALUE OF THETA = 
INCREMENT OF THETA = 


100 MAGNITUDE OPTION = STD 
. OOOOO0E+00 
s2oueOOR-01 
ten UTERATOR CO. “FICIENTS bC0O),b( 1)... bCE) ARE: 
© 2900E+00 nme rac 2 lei 0) 


- 8811E+00 maZ00E+00 


ibe PeNewINATOR COEFFICIENTS cC0),¢c(1).2360CN) ARE: 


| 


. LOOOE+01 . OOOOE+00 moe WZ E700 . 3234E+00 
OUTPUT DATA FOR SYSTEM # 1 
THETA MAGNITUDE PHASE 
(RADIANS) (DEGREGS: 
. OOOOODE+00 . 999940E+00 . OOOOODE+00 
WZ OOO00E-01 . 999940E+00 =~ LOS OUR On 
fee 2005-01 £9999 39E 500 aes OO 29 S01 
. 841600E-01 mg 29 39000 moO 0432650 | 
. 112240E+00 79999 39R.,00 aa (9 / 7S On 
. 140300E+00 9999 36E700 as 77 Ook a0. 
. 168360E+00 999957 EO0 L177 SOEL00 
. 196420E+00 799995 0E00 oe OOO Se O0 
. 224480E+00 2 9999S5E700 ee to O> 7) 9Ee.00 
722 0E+00 . 999934E+00 ol / 947 SR eO0 
- 280600E+00 Bel elehe 6) Sec d0lU =, 200768E+400 
. 308660E+00 ,999932E400 ~e222) loEmuO 
yeeo77 OEt00 9999 30Es00 ~. 244780E+00 
. 364780E+00 . 999928E 500 = 207/018E,00 
. 392840E+00 7 9999ZCE 00 oe SEN O10, 
- 420900E+00 9999 25E00 =. 01 5332E.,00 
. 448960E+00 999921E 00 -. 340364E+00 
. 477020E+00 D299 9TBETOU = 300305E7,00 


- SUSUS0E00 
. 933140E+00 
, SOlZOGE-U0 
67 200E 00 
7017 g20n 700 
. 645380E+00 
. 673440E+00 
2/015 008,00 
2/7 29560E200 
75/60Z20E200 
. 785680E+00 
. 813740E+00 
. 841800E+00 
. 869860E+00 
Baie) SVAN ON8, 
- 925 950E7T00 
. 954040E+00 
- 9621005700 
- ICLOTGESO1 
« LO0622E01 
2IOGG26E.01 
. 109434E+01 
.112240E+01 
. 115046E+01 
pl? So ZEO 1 
wiZOG>6Ea0 1 
. 123464E+01 
22627 OEE0 1 
, 2907 6E+01 
odes Zia 0) I 
-os056n701 
.137494E+01 
. 140300E+01 
Pres TOCET OT 
a AS AT Oa 
.148718E+01 
p hohs 4 an Ol 
-2 54 5008701 
#157 1Oon.01 
5 US eee sOUl 
. 162748E+01 
- 165554E+01 
PLoS SC OE. Og 
shy alleen 2c 40) 
eye Ade AO 
e767 7ou.0d 
1795648501 
~ 18239080 1 
eS De ioest 0 1 
So OO 24-01 
SIOSUSE 4-01 
seo oO ert Od 
. 196420E+01 
WAZA ead0ul 
PeOZ03ZE 01 
. 204638E+01 


» 9999 1400 
og ee 0) 
STS NOs 9).0) 
. 999902000 
Bele retell 1@) 
-ISDCITETOO 
. 999884E+00 
999877 a0 
. 999868E+00 
99965 955,00 
. 999848E+00 
IIDC SSEROO 
2997 8c2UE OC 
- 9998038200 
. 9997 E3E500 
19997 59-00 
9997 SEaOO 
. 9996097 E00 
Rae we los) ota 010, 
- 99960er +00 
- 999>s SE 00 
-999454E+00 
: 999642500 
. 99D TICE TOO 
Belehercleo ean 00, 
; 998631500 
99627 ers 00 
.997441E+00 
» 996025E-00 
s 992995 ,00 
. 984532E+00 
- 942860E+00 
- 38203 2E On 
, 999922E2,00 
. 984110E+00 
- JI 2860E FOU 
99597 0H Eo 
99/7 4138s 00 
, 99 6 O70 
. 99367 1E00 
. 99697 SE Ue 
99° [GOR ee 
- 999340E500 
- 999452E700 
. 999537Es0e 
. 99960300 
- 999655E 700 
» 99969 7E0G 
» 99973 2na0e 
»999701E2 00 
999 /65E200 
, 999805E5,00 
. 9996 27E eG 
- 999837 Hae 
- 999850E+00 
. 9998ol1EEO9 


-29 525 000 
4213278200 
-450665E+00 
- 481406E+00 
olay 2 100 
-547805E+00 
wo 026 (5b t00 
- O22183E500 
S665 03757-00 
. 706764E+00 
7237 OoEro0 
. 804590E+00 
. 859746E+00 
- 919964E7+00 
 SBC1L0IETO0 
, JOS 922 ou 
- 114061E+01 
. L231 59E+ OL 
. §35540E7,0 1 
. L45369Re0n 
- loJeo>E ron 
»lj7ou7 oErO1 
. 194344E+01 
- 217625 E01 
. 246560E+01 
i 2O0D DEO. 
Hoes Vda sO 
. 400968E+01 
, S020 5ens01 
0” 29601 
025 ORR 2 
» 194058E+02 
. 67 7 So4E 02 


PAGO SIs a0) 
s TO2Z59SE G02 
; 68/1568 .01 
22 LOLI 201 
»412999E+01 
. 344004E+01 
. 294544E+01 
5 )7) c 10K 16 or (0) 1h 
- 228225001 
. 204850E+01 
185627501 
+ LOIS ZUROd 
- 155807E+01 
-143977ErOl 
Slee los 718 210/11 
~124549E+01 
- 116454E+01 
el OS VOT Ea oa 
. 102646E+01 
. 966947E+00 
29125665700 
~OOZ621EOW 
eo loo Z22Ee9 


OZ 


.207644E+01 
. 210450E+01 
fee OCrE +O) 
melo Oe 2-01 
eciSsoosk+0 | 
.221674E+01 
. 224480E+01 
Poe OorO | 
e222 009ZE70 1 
eee eons] 
2 255)/7046E+01 
pZo5 5 108501 
P2eS16E+01 
.244122E+01 
w2eo7 ZOE 1 
» 249734E+01 
wee co408 +0 1 
eS oor) 
256 15 ZE701 
720075087,0 1 
. 263764E+01 
me Ooo) Ob 0 1 
weed. 7 Obs 0 | 
ee rea OL 
. 274988E+01 
me 77 IGE | 


, 9998 JAE SOO 
» I9F5S0E 00 
- 999887 7E00 
- 999894E+00 
e229 970Ch 00 
. 299 90EE.00 
99971 Teo 
79999 1s EO 
22999198 00 
. 9999 23u 7,00 
. 9999268 2,00 
9999 29E 00 
. 999932E+OD 
» I9DISGETOO 
. 999 SOE LOO 
. 9999S9 00 
9929 C0E +00 
. 999942E+00 
. 999944E+00 
~9 99945 E700 
. 999947E+00 
.999948E+00 
799992300) 
: 9999S0Es00 
999905 LE rae 
, 99995 2200 


J J7S7G56 700 
. /34003E+00 
. 696833E+00 
. 661978E+00 
76296 ve 00 
Beye lic\4 Ss) 0g OL8. 
266958 E500 
Sai s00 
. 514669E+00 
- 489423E+00 
, G052 778500 
-442138E+00 
.419907E+00 
. 398514E+00 
nov? O7 OF LOU 
moo 3700 
To 55029800 
5199 05Ec00 
- 3017 168200 
. 284016E+00 
. 266763E+00 
» 249920E+00 
. 233452E+00 
acl SZ26E500 
220 15078700 
- 185980E+00 


END OF RUN, SYSTEM # 1 


D. MAGNITUDE, PHASE OF THE AUGMENTED IIR NOTCH FILTER: P= 


INPUT DATA FOR SYSTEM # 1 


INPUT DATA SOURCIFILE: digfreq. tst 


DEGREE OF NUMERA™ OR = 4 
DEGREE OF DENOMINATOR = 4 
dsorce =F 
NUMEER OF FREQUENCY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = . 000000E+00 
INCREMENT OF THETA = .280600E-01 
THE NUMERATOR COEFFICIENTS b(0),b(1)...bCL) ARE: 
. 9900E+00 -. 3299E-02 . 1861E-01 .6141E+00 -. 8625E+00 
THE DENOMINATOR COEFFICIENTS c(0),c(1)...c(€N) ARE: 
. LOOOE+01 . OOO0E+00 . OOOOE+00 . 6108E+00 - 


. BSSIETOO 


oe 


Re 


HEA 
(RADIANS) 


. OOOOO0E+00 
-2o0G0Gr-O1 
sDoleeOr = O01 
.841800E-01 
. 112240E+00 
. 140300E+00 
. 168360E+00 
. 196420E+00 
. 224480E+00 
» 252540E+00 
. 280600E+00 
. 308660E+00 
, 3507205100 
. 364780E+00 
. 392840E+00 
. 420900E+00 
. 448960E+00 
-477020E+00 
TOSCO OUE OO 
- 933140E+00 
; DolZO0E00 
. 589260E+00 
toll 9 72. 0)3.2 016) 
. 645380E+00 
. 673440E100 
. /O1S0CE+00 
w/29560ES 00 
so TO cUnTOU 
w/B>6GUET00 
.» 813740E+00 
. §41800E+00 
. 869860E+00 
Beep ih P10) a 1018, 
92595 0ETO0 
. 954040E+00 
27012 100E 00 
s1O1OT6E+0 1 
pile eZ ZEr0 1 
mi ooZGrs-0 1 
. 109434E+01 
»112240E+01 
. 115046E+01 
ay Oo Zod 
 IEZOSCS OE Om 
. 123464E+01 


MAGNITUDE 


. 999940E+00 
. 999940E+00 
. 999940E+00 
. 999940E+00 
»29I9D ISI EOO 
- 9999 36000 
- 9999 SC EEOo 
. 999937 ETOS 
. 9999 S6E00 
. 999934E+00 
©9999 S3En00 
29799 3 2a 
. JIIISCETOS 
Gee he 70) aA 0h0. 
9999 2ZSE TOO 
» 99992 SEs OO 
99722 OLEOO 
Re SL A0i6 
ooo oO Lees OO 
29975 FOnT00 
99920 600 
Me, ele) Olin 018, 
F elelehey slots 306) 
» 9922915. +00 
. 999884E+00 
99957 6n4,00 
. 999868E+00 
2 2 220 5CE 00 
1999 347ET00 
R cle iets Ss) dad 016, 
ee iio (210) Sa AON 6, 
Be 2!) (o10 S)o5 1010, 
S999 7 63 EsO0 
“Posi SoEtOo 
9997 S150 
Oly EEO 
Se) SS Sdn 1010) 
; 99960 5E100 
5 99953 7200 
- 999453E+00 
,999Sa5n 700 
2999 193E.00 
727627 65E..00 
- 9986G1E00 
»7962 lsEs UC 


OUTPUT DATA FOR SYSTEM # 1 


PHASE 
(DEGREES ) 


. OOOOO0E+00 


94 


, 194072E- 01 
wooc4Zon-u ll 
=553013n-01 
7 SOC 22Ee01 
577 O0souael 
FIL SOZECe 
. 138056E+00 
. 158606E+00 
,179500E 00 
. 200794E+00 
» 222545E+00 
» 244805E+00 
. 26/638E+00 
27 LIS6hreo 
, SIS 351ET 00 
. 340386E+00 
, 360521 E00 
~ooo27 UE00 
»421345E+00 
-450679E+00 
-481422E+00 
5137 34h 00 
»54/7815E+00 
- JOJ005 E00 
, O22a7 E00 
. 663044E+00 
s OG 7SE00 
| 7537 96E.00 
. 804595E+00 
, O59 Jobo 
.919967E+00 
. J86NOZEROG 
2 105 921E Ol 
. 114062E+01 
Vitale kop bg 
.133548E+01 
. 145388E+01 
159094EF01 
ule ll 7 SiO ll 
. 194343E+01 
217 625Ea50 Il 
. 2465568701 
5203257 E40 1 
So o2Os0h+OL 


mlZe2 20E+01 
ey OE +O | 
mele oo7 E01 
. 134688E+01 
milo 4 EO!) 
- 140300E+01 
. 143106E+01 
pies 12E+01 
. 148718E+01 
ese 26h 0 | 
p45 S0E+01 
fey Loon+O1 
mile 74 2E +0 | 
. 162748E+01 
wlees 5460 | 
mlocso0H+01 
227 P166R+-0 1 
slyoo7ZErOl 
miO7 7 6E0 1 
mio 6 Gr 01 
J tend SPT OP Oe OM 
ee Or 0 
» L6c002ZE+01 
. 190808E+01 
pleco ler +o 
et 76e20E+01 
pee? 26r0 1 
powewo7 Lao | 
. 204838E+01 
. 207644E+0]1 
. 210450501 
» 2132560 1 
P2160e2ir C1 
eoleseouro | 
.221674E+01 
. 224480E+01 
pees Zoors UL 
eZo009ZE+01 
eee ov onsO 1 
yeoe /OCE+0 1 
aoe 5 LOR 1 
-2o8o160E+01 
.244122E+01 
a2a0925E+01 
-249734E+01 
.252540E+01 
.255346E+01 
ee ere Bo) | 
eeu IDOE 1 
. 263764E+01 
ge ceo OL-OiL 
Pogo 27 OE0 1 
5 1) Me eB dns OM 
Pe7e3e ert 1 
e277 7 94E+01 


- 997440E+00 
; IO O 245-00 
199299 GEE O0 
- 984530E+00 
. 942856E+00 
eves s 7120 
. JS99ZCEFOO 
270411 1E+00 
$992 007 E00 
Ree 3) Sie Opn 80) 
IIT O13E OO 
S99 e196E E00 
- 72007 LE00 
J920979E 100 
Bele allies os )0 
9993405200 
9994525400 
ooo? bOw 
. 99760 3E+00 
2 IPICOS ao 
IOS EEO 
s IS oF 32E 00 
ieee OlEToO 
-999/755E+00 
29276 0on00 
2999 622E00 
wg goo7 bo” 
7227600200 
2999COLEOO 
799987 FETO0 
yee sOreo 0 
. 999 6o7 hs-00 
- 999694500 
- IP29COE TOO 
799 9906E+00 
Sele ee eas O18) 
fe) Sle Ske tan N6. 
eee eee ES Oi 810) 
« 99922 EO 
so 9 Zon 00 
29999298200 
Res) el el leads 0) 0. 
; 9999348700 
. 9999360200 
e999 9352700 
. 999940E+00 
799994 2500 
. 999944E+00 
. 999945E+00 
. 99994678400 
. 999948E+00 
» 999949E00 
: 9999502700 
99995 Et Ow 
JI 909 SZELOe 


° 


400966E+01 
502849E+01 
O72 3 6b Oi 
100258E+02 
194057E+02 
S77 so35 8702 


20001 ESOZ 
»LlO2Z2594E+02 
Corts lEOd 
, SISTOSESO1 
~olZ595E+01 
. 344001E+01 
. 294542E+01 
487 Sloss Ol 
yeeoecous Ol 
. 204849E+01 
1 lesoZ6naol 
5 kisnehsy Ikeyoiar O) Il 
PIS >OUGE+-O1 
mle eo) Oe to 
7 ieee > Ol go)L 
.124548E+01 
- 116454E+01 
OL Ion +O 1 
. 102646E+01 
slots esp oa5 O10, 
- 912560E+00 
Poocolot.e0 
sels Silica 010, 
777 7 EEOC 
. /34001E+00 
ogee ZE+TU0 
~OOl77 7 E00 
629 1SZET00 
~peoeo TE TO0 
0c o 20 L700 
-541134E+00 
.514669E+00 
- 489424E+00 
THO 27 9EtOU 
.442136E+00 
2 abo 910E O00 
Seon Aue 
Bee Is asda Ole, 
oo) Jeo OU 
wasoOS2 E400 
pe Seca Oe 
OO f7 reno 
. 284017E+00 
~200705n 1.00 
. 249920E+00 
. 233450E+00 
52 IES Vardi gl 818 
5 ADS Boas Ole 
lis eke re pda5 00 


JIS) 


E. MAGNITUDE, PHASE OF THE AUGMENTED IIR NOTCH FILTER: P=3 


INPUT DATA FOR SYSTEM # 1 


INPUT DATA SOURCEFILE: digfreq.tst 
DEGREE OF NUMERATOR = 5 

DEGREE OF DENOMINATOR = 5 

dsorce = F 
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NUMBER OF FREQUENCY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = .0Q0QQ00E+00 

INCKEMENT OF THETA 2” = 2806G00E=01 
THE NUMERATOR COEFRPICIENTSEbCO) (DG)... =) (Gneee in 

JI SVOETO 0 -. 3299202 wIS6LE=01 .9374E-02 =. 6610400 
-. 6047E+00 
THE DENOMINATOR COBFFICIENTS =eU@)) ye C1) ee cree: 

SOCORRO} . QOOOE+00 . OOOOE+00 . OOOOE+00 =, 65245706 
S57 6/7 a0 

OUTPUT DATA FOR SYSTEM # 1 
THETA MAGNITUDE PHASE 

(RADIANS) ( DEGREES ) 

. QOOOQ0E+00 - 999937500 . QOOO000E+00 

| 2o0500R-01 JI9DOSSEsO0 -. 194494E-01 

wo CmZOUE S01 - 999938 OO = 9069 AUB On 

. 841800E-01 - 999938E2,00 -.584018E-01 

- 112240E+00 . 9999SSE 00 =. /80285 011 

. 140300E+00 meee iors) an070 =. 9/815 70201 

Sl etee oO dab) 8) 999937 RE LCe = 117817 E00 

. 196420E+00 - 999936E500 = I3S06an 00 

. 224480E+00 - 999935nE00 =, Logor0r.00 

20 2040E00 / 999FIS4E LOO =, 179497 Eee 


. 280600E+00 
. 308660E+00 
ma20/Z0E+00 
. 364780E+00 
poo 2o4O0E+TOO 
- 420900E+00 
. 448960E+00 
.477020E+00 
eS 050S0E+00 
peoe 140E+00 
wo 200E+00 
mpG2200E00 
ome oZ0E+00 
- 645380E+00 
- 673440E+00 
2/7 O1500E7T00 
#7 29500B+00 
mie OZ 0ET00 
He OCORLUD 
. 813740E+00 
. 841800E+00 
. 869860E+00 
pod) J20E+00 
po2599805400 
. 9540409E+00 
62 100E400 
nO TO16GE-01 
mero 36228401 
~lUo6ZGE+01 
. 109434E+01 
Pere 26 OE 0 | 
| 15062701 
mile os 225,01 
wi206565+0 | 
. 123464E+01 
mo 2 Obs Ol 
lz 907 OFT 01 
Pole 2b 01 
. 134688E+01 
2 374G05E+01 
. 140300E+01 
ple S1OCE+01 
poe 1 2E+0) 
- 148718E+01 
eso rt O 1 
7i54330E+01 
ils FObSiopop soul 
love 2h+0 1 
ml62746E+01 
. 165554E+01 
. 168360E+01 
floor +0 | 
ee 397 Ze) 
ml o7 7 OE+01 
sly o5o4k+01 
ele2o70ETO! 


, 99993 3R400 
99993 TES 
- 999929 
3999927E 00 
eI IIS 2S Ee 
» 99992 SE Oe 
» 999S20E+00 
599997 E00 
. 9999 lak 00 
. 9999 10E TOG 
, 9IFFUCE 00 
99990 2ZEOw 
- 999896200 
~I9SISI LESS 
. 999884E+00 
- 9998 76E700 
, 99D S6GE 00 
999 65SE200 
99984 /ELOU 
pele eho pe ied 3010, 
7 999S2Z0E 00 
, I9ISOSE AIG 
-999783E £00 
/ 999 /S5SEr00 
. 99973 1EROO 
_ 999697 E200 
999655500 
"IPP COSE TOO 
IIIS Sy eo 
79945 SE 00 
. 9993422400 
999 Toei 
»996965E500 
“99665 TETOO 
29902 WgEroU 
. 997440E+00 
. 996024E+00 
~99Z9945E+O0 
. 984530E+00 
- 942856E+00 
sooo 35k =01 
pg 399Z45E+00 
. 984110E+00 
Hee2ZO007E TOO 
29959705400 
; 9974613E+00 
Boo OLI6E 00 
799867 LETOO 
SI9Oo7IETOO 
nooo lee O00 
. 999340E+00 
7 999452E7,00 
nodD OOO 
eV gOOUSE TOO 
Heo] 5E7 00 
poo geo) E00 


- 200787500 
- 2225305000 
. 244798E+00 
+26/6355200 
29 bI27EaOe 
. 315344E+00 
3403/7 5E500 
Spoleto 2)3 2210) 9, 
PO93Z6Sr 00 
.421337E+00 
st 06/7 1E +00 
-481411E+00 
fo 13727600 
-24/810E+00 
,205876E700 
- 622189E+00 
P6635 05eE7.00 
. 706766E+00 
. 753794E+00 
- 804590E+00 
.859744E+00 
ao 1SI6IEs00 
77 SOU OE. 00 
51059 2b Or 
»114061EF0m 
a2 2199E 01 
mis 3547 EroR 
. 145388E+01 
eo 9094 Es Oa 
PITS IT SETO 
5 ESTES NiO) Th 
m2 1/7 625E701 
. 246558E+01 
we OSD oy EaOl 
moc oo Or a. I 
. 400966E+01 
. 502848E+01 
Ov IZ SGETO! 
SIOOZS58E7,02 
ge 0so7EaO02 
Of eS O2 


P2000 TIE OZ 
lOZ5946Es02 
moo? 15 1B LO 
-516108E+01 
eG 12995E+01 
. 344001E+01 
eons OIL 
5, 5)) SOoldnd eal 
2262235401 
. 204848E+01 
7 o> OZOn+01 
moo leet Od 
- 155605E701 
militeo 97 SE +01 
mils2 65 0E7,01 


oy) 


~ISSLSCE+Cu 
_- L6souZznst or 
J 90S0SEa0" 
. 193614E+01 
; Wea Zons oa 
1992 26E701 
220202 2E 00 
-204636E401 
.207644E+01 
7204 S0RG 
eZ S25 or ave 
~216062h 501 
5 Pere leved dar Oil 
222 67 GeO 
. 224480E+01 
wee 2O0nrUL 
2 500C2E-01 
PSP oys) Blas OL 
725570GE201 
5 sehe SOI e (0) Il 
- 241316E701 
» 244122E+01 
- 246928E+01 
. 249734E+01 
e252 5008 0d 
72> Soar +Or 
eo 2 
220095 Ch704 
22057 OF EOL 
e200 57 0be0 1 
720957 Oc ou 
72/2 1N62E Un 
27 E9SOR40 I 
27 oO 


7999732500 
-99976ar oo 
- 9997 35E-00 
9 I9805ECOe 
99982 2Emu0 
. 999637 ueCe 
: 99SSS50E+00 
+999661Es.00 
7999 67 ee Oe 
- 999 SCOEECE 
- 999S87ESO0 
. 999894E+00 
- 999900n ade 
> 999 90GE 00 
~ 9999 TIETO0 
-IDIGISESUG 
- I9DDISETCS 
- 9999275 BOW 
- 9999 26R Ue 
- 9999 25E 00 
/ 99993 2E700 
79999348700 
- 999936E 00 
,999936En00 
~9999L0ES eo 
-999942E+00 
- 999943E=00 
- 999945E200 
. 999946E+00 
: 9999SSE Te 
,999949R-00 
, 999 950En00) 
- 99995 1E208 
- 999952500 


. 124548E+01 
- L16453E+01 
1091 96r ot 
. 102645E+01 
7 906095/7E ie 
-912553E—09 
/ 66260GE00 
761605068700 
2/73//0E 0G 
2 /335993E200 
7696525E200 
. 661968E+00 
SO2Z9177E +00 
7>98225E 200 
7 508932E-0c 
~541 1258406 
. 514664E+00 
-489416E+00 
-465272E+00 
- 4421298406 
.419903E+00 
; 3598506E-00 
.3//6/2E500 
795/93 1E neo 
, 53803 1ES00 
oS LOYT Eo 
. 90172TESO0 
. 284018E+00 
72660771 EeOe 
- 249930E 700 
. 233465E+00 
.217347E 560 
»201533E500 
- L8600ZE—00 


END OF RUNG Sars ire 1 
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APPENDIX B. PROGRAM AND DATA OF DIFFERENCE EQUATION 


A. 


AAQAAAAAANAANGAAAIANAIANAIAINANAIANAIAIAIAIIAIMIAIAIAIMIAIIAANIAIIAIAIaAGY 


Pan! 
¢ 


PROGRAM FOR DIFFERENCE EQUATION CALCULATION 
Ee OO VERSION: 2/03/88 


PURPOSE: THIS PROGRAM COMPUTES THE ITERATIVE SOLUTION TO A 
LINEAR, TIME-INVARIANT (LTI) DIFFERENCE EQUATION. 
THE DIFFERENCE EQUATION MUST BE IN THE FORM: 
y(ns) = a(i)*“y(ns-1) +... + a(N)*y(ns-N) + 

b(0)*x(ns) + b(1)*x(ns-1) +... + bCL)*x(ns-L). 

THE PROGRAM CONSISTS OF A MAIN PROGRAM AND TWO 
SUBROUTINES. SUBROUTINE diffeq IS CALLED BY THE MAIN 
PROGRAM TO ITERATIVELY SOLVE THE DIFFERENCE EQUATIONS AND 
SUBROUTINE xgen ALLOWS THE USER THE OPTION OF GENERATING 
THE INPUT SEQUENCE x() BY WRITING THE APPROPRIATE 
EQUATIONS. IF THE USER ELECTS TO GENERATE THE SEQUENCE 
x() BY SING xgen THEN THE PROGRAM MUST BE COMPILED 
AGAIN BEFORE EXECUTION. THE USER HAS THE OPTION OF 
SELECTING ONE OF TWO OPERATING MODES: BATCH OR TEST. 
IN BATC! MODE THE AMOUNT OF INTERFACE WITH THE USER 
IS MINIMIZED AND IT IS ASSUMED THAT THE INPUT DATA 
HAS BEEN STORED IN THE DEFAULT FILE 'DIFFEQ.IN'. IN 
TEST MODE THE USER IS PROMPTED FOR THE NAME OF THE 
INPUT FILE OR HAS TEE OPTION OF PERFORMING A TEST RUN 
USING THE INPUT DATA STORED IN THE FILE 'DIFFEQ.TST'. 
IT IS RECOMMENDED THAT FIRST-TIME USERS SELECT THE 
TEST MODE AND MAKE A TRIAL RUN WITH THE PRESTORED 
INPUT DATA. THE TEST MODE ECHOES PORTIONS OF THE 
INPUT DATA ONTO THE MONITOR TO ALLOW VERIFICATION OF 
ITS ACCURACY. BOTH BATCH AND TEST MODES ALLOW THE 
USER TO SOLVE UP TO FOUR DIFFERENCE EQUATIONS IN A 
SINGLE PROGRAM EXECUTION. THE OUTPUT OF THE PROGRAM 
'DIFFEQ. FOR' IS STORED IN THE ARRAY y(). THE OUTPUT IS 
STORED IN TABULAR FORM IN THE OUTPUT FILE 'DIFFEQ. FOR' 
AND IN A FORM SUITABLE FOR PLOTTING IN THE FILE 
"DIFFEQ. DAT’. 


ee aa eR ee [Peet eee ed ee ey ae ew Peel eel ee eee, lo ale Santectwaloatoa! wWoulscloataatoctout, ’ ee aa San wl, fawewewewantantan! 
Soe ei bari bir Cacion facie | oketcarl ok low bo ir ok oir tok tok Conroe bora ity INPUT sesovesevevesovesesese se dese sevesesesesehete tated stele eked: 


THIS PROGRAM ASSUMES THAT EACH DIFFERENCE EQUATION IS IN THE 
CANONICAL FORM: 


hems) = ac l)<y(ns-1) + .:. + a(Njeytns-N) + 
b(0)*x(ns) + b(1)*x(ns-1) +... + bCL)*x(ns-L) 
A NON-NEGATIVE INTEGER, THE NUMBER OF INPUT DELAYS. 


a 
Hol 


A NON-NEGATIVE INTEGER, THE NUMBER OF OUTPUT DELAYS. 
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a(1)...aCN) = REAL COEFFICIENTS OF THE GUPPUT TERMS: 
b(0)...b(L) = REAL COEFFICIENTS OF THE, INPU tERMce 


THE INPUT PARAMETERS SHOULD BE STORED IN A FILE NAMED 
'DIFFEQ.IN'. ALL OF THE READ STATEMENTS USED BY THIS PROGRAM 
REQUIRE FORMATTED INPUT. PARTICULAR ATTENTION SHOULD BE PAID 
TO THE FORMATS; ESPECIALLY THE USE OF THE SOER ie. o ni ne 
DENOTE 'REAL’ NUMBERS. THE INPUT PARAMETERS REQUIRED BY THE 
PROGRAM ARE LISTED BELOW. 


NAME TYPE RANGE (ARRAYS) RESTRICTIONS 
numsys PINTER 1 <= numsys <= 4 
L INTEGER O <= L <= 128 

N INTEGER O <= N <= 128 
nstop INTEGER QO <= nstop <= 999 
xsorce CHARACTER ORS 
b(k) REAL OLS cose ane O <= L <= 128 
a(k REAL Laos 3) eeeeN O <= N <= 128 
y(k) REAL =) eee | 1 <= N <= 128 
x(ns ) REAL OL ly -aeomstor QO <= nstop <= 999 
WHERE: 


numsys = THE NUMBER OF SYSTEMS TO BE EVALUATED. 
THIS INTEGER VALUE MUST OCCUR AT THE TOP OF THE INPUT 
FILE. IT DELINEATES THE NUMBER OF SYSTEMS TO BE READ BY 
THE PROGRAM AND ANALYZED. FOR EACH SYSTEM (1...numsys) 
THE PARAMETERS BELOW MUST APPEAR IN THE INPUT FILE. 


L = AN INTEGER VALUE THAT SPECIFIES THE MAXIMUM NUMBER OF DELAYS 
_IN THE INPUT SEQUENCE. 


iz 
1 


AN INTEGER VALUE THAT SPECIFIES THE MAXIMUM NUMBER OF DELAYS 
IN THE OUTPUT SEOUL NCE 


nstop = AN INTEGER VALUE THAT SPECIFIES THE LARGEST TIME INDEX 
(ns) FOR WHICH THE DIFFERENCE EQUATION IS TO BE SOLVED. 


xsorce = A CHARACTER VALUE OF 'F' OR 'S’ DENOTING WHETHER THE 
INPUT SEQUENCE x() IS TO BE READ FROM THE INPUT FILE (F) 
OR TO BE GENERATED (S) USING THE SUBROUTINE xgen. THIS 
LATER OPTION IS ATTRACTIVE WHEN nstop IS A LARGE NUMBER 
AND THE INPUT SEQUENCE x() CAN BE READILY DESCRIBED BY AN 
ANALYTICAL EXPRESSION. IF xsorce = 'S'’ THE USER MUST 
PROVIDE THE APPROPRIATE FORTRAN STATEMENTS IN THE SPACE 
PROVIDED IN SUBROUTINE xgen AND THE PROGRAM MUST BE 
RECOMPILED BEFORE EXECUTION. 


b(k) = REAL COEFFICIENTS OF THE INPUT SEQUENCE x(ns-k) IN THE 
ORDER: b(0), b(1), ..., b(L). 


a(k) = REAL COEFFICIENTS OF THE OUTPUT SEQUENCE y(ns-k) IN THE 


ORDER: a(1), a(2), ..., a(N). IF N = 0 THEN THE EQUATION 
IS NON-RECURSIVE AND NO a(k) COEFFICIENTS SHOULD BE IN 
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y(k) = 


SUE SINEUE FILE. 


THE INITIAL CONDITIONS FOR THE OUTPUT SEQUENCE IF THE 
VEE PERENCESEQUATION IS RECURSIVE, 1:E., N > 0. THIS 
PROGRAM CALCULATES THE SOLUTION TO THE DIFFERENCE 
EQUATION FROM ns = 0 TO ns = nstop THEREFORE THE INITIAL 
CONDITIONS y(-N) TO y(-1) MUST BE PROVIDED IN THE INPUT 
PE Ne tnemenhOnR: yC=-N)seyC-Ntl), .sageiyG-1). IF N = 0 
THEN THE EQUATION IS NON-RECURSIVE AND NO INITIAL 
CONDITIONS SHOULD BE GIVEN IN THE INPUT FILE. 


x(ns) = THE INPUT SEQUENCE. IF xsorce = 'F’ THEN THE INPUT 


SEQUENCE x(0), ..25)x(nstop) MUST BE PROVIDED BY THE USER 
IN THE INPUT FILE. IF xsorce = 'S' THEN THE USER HAS 
ELECTED TO GENERATE THE INPUT SEQUENCE BY PROVIDING THE 
APPROPRIATE FORTRAN STATEMENTS IN THE SUBROUTINE xgen. 





NOTE: THE INPUT FORMAT STATEMENTS OCCUR IN THE MAIN PROGRAM 
FOLLOWING THE CAPTION:  xdsdededsdece INPUT FORMAT =o wiseevedededede | 
THE FORM OF THE INPUT DATA FILE IS: 
LINE # ENGR TES FORMAT 
1 numsys a 
Zz L,N,nstop,xsorce deoeee aor 28 leona Oar. 
NEXT NB LINES Gian ee Oy eee ols 6£10.0 
Ext NA LINES AG es NN 6£10.0 
Bent NY LINES VC a hae ee gtk 6£10.0 
NEXT NX LINES x(ns), ns= 0,...,nstop 6£10.0 
WHERE: NB = 1 + (L/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 
NA = 0 IF N = O OR 
NA — 1 + ((N-1)/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 
NY = 0 IF N = 0 OR 
NY = 1 + ((N-1)/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 
NX = 0 IF xsorce = 'S’ OR 
NX = 1 + (nstop/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 
MEGxsorce = sce 
*NOTE: FOR numsys > 1 THE FORMAT OF LINES 2... IS REPEATED. 


THE FORM. T £10.0 USED FOR INPUT DATA PERMITS THE DECIMAL 
POINT TO BE PLACED ANYWHERE IN THE FIELD OF 10 COLUMNS 
AND ALSO ALLOWS THE EXPONENTIAL FORMAT TO BE USED (EG. 
3146.2 = 3.1462E+03). 


Sevedededesvedededvesevesededededededesedevededededeslededeses OUTPUT sksedesecesteskscsestescsescsestededesevevescckedese cle se deck deve 


THE INPUT DATA AS WELL AS THE OUTPUT DATA ARE STORED IN TABULAR 
FORM IN THE FILE 'DIFFEQ.OUT'. ADDITIONALLY, THE INPUT SEQUENCE 
AND THE OUTPUT SEQUENCE ARE WRITTEN INTO THE FILE 'DIFFEQ.DAT' TO 
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FACILITATE PLOTTING BY A SEPARATE, USER SUPPLIED PROGRAM. THE 
FORMAT OF THE DATA IN 'DIFFEQ. DAT’ IS: e12.6, 2X, e12.6. THE FIRST 
ENTRY CORRESPONDS TO THE ORDINATE VALUE, AND THE SECOND ENTRY, THE 
ABSCISSA VALUE. ADDITIONAL HEADER INFORMATION IS WRITTEN INTO 
'DIFFEQ. DAT' TO ALLOW FOR CONTROL AND LABELING OF EACH PLOT. 


mer ae orice icy totrige ise ye aceite sine ice pris eee ses rises EXAMPLE seseseosscskesesesedeseseseocsesteesesesesesesdkevescdcoesese 


THE INPUT PARAMETERS FOR THE SYSTEM DESCRIBED BELOW ARE STORED IN 
THE SAMPLE INPUT FILE 'DIFFEQ.TST' AND CAN BE USED FOR A TRIAL 
RUN IN TEST MODE. 
DIFFERENCE EQUATION: 
y(ns) = 1.2*y(ns-1) + 1.5*x(ns) 
GOAL: TO OBTAIN THE SOLUTION TO THIS DIFFERENCE EQUATION FOR 
ns = 0 TO ns = 10, GIVEN: x(0)...x€10) = 100.0 AND 
THE INITIAL CONDITION y(-1) = 25.0. 
THE INPUE Files: 


1 


000 001 010 13 

Lo 

Le? 

Zon 

100. 0 100.0 100. 0 1CO..9 ILO} 0, 10 he O20 
POOR 100-0 100. 0 100.0 10020 


THE RESULTING OUTPUT FILE 'DIFFEQ. OUT' IS: 
INPUT DATA FOR PROBLEM # 1 
PROBLEM # 1 INPUT DATA SOURCEFILE: DIFFEQ. TST 
THE NUMBER OF INPUT DELAYS: L= 0 
THE NUMBER OF OUTPUT DELAYS: N= 1 
THE VALUES ORES tes ls va nC 
THE COEFFICIENTS BUC) 9 bG))eeeEGny ARE: 


. LOOCOUCR Re 


THE COEFFICIENTS a(1), ..., a(N) ARE: 


. 12ZOO00E+01 
OUTPUT DATA FOR PROBLEM # 1 
ns xCns )} y(ns) 
-1 . 0000. JE+00 . 250000E+02 


102 


0 . LOOQODE+03 . LSCOCCE:03 
Hi - 1000CCE+03 2 66000E+03 
2 . LOOOOOE+03 sD 6920003 
3 . LOOOQOE+03 .857040E+03 
+ . LOOQOOQE+03 .117845E+04 
5 mEOCCUCES03 . 156414E+04 
6 . LOCOCOES03 . 202697E+04 
7 mEOCOOUETOS »2oO250E+04 
8 . LOOOOOE+03 . 324883E+04 
9 Zo VOOCETOS . 404860E+04 
10 mECOCOOETO3 . 500832E+04 


eee END OF PROBLEM # 1 ---------------- 


FOR ILLUSTRATIVE PURPOSES THE INPUT SEQUENCE x() COULD HAVE BEEN 
GENERATED BY SPECIFYING xsorce = 'S' AND WRITING THE APPROPRIATE 
FORTRAN STATEMENTS INTO SUBROUTINE xgen. THE STATEMENTS THAT 
COULD BE USED TO ACCOMPLISH THIS ARE WRITTEN INTO THE SUBROUTINE 
BUT ARE ‘COMMENTED OUT'. 


IAMIAAIAIA’AIAIAINIIMAIIAIIAIIMlI|IaIIIa|Ia 


5 
foun eelse 20 nea e's we'wo en wale suse sales clones eens a'eae att - one OY Veh ak Veletal . - ae'ee se oe ee ete we ers 
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character infile*12, mode*l, xsorce*l 
meet 2S), ClOL1286) ~y(-126: 850), x(-126: 850), a2 


fee ee PT eUSEK FOR MODE: BATCH OR (TEST. 


write(*,1115) 
read(*,1117) mode 
if((mode. eq. 'Y'). or. (mode. eq. 'y' )) then 
mode =  Y° 
write(*,1118) 
read(*,1119) infile 
else 
infile = ‘DIFFEQ. IN’ 
endif 


C UNIT=1 DEFINED AS INPUT FILE. UNITS=2,3 DEFINED AS OUTPUT FILES. 


open(unit=1,file=infile,status='old' ,iostat=ierr,err=999) 
open(unit=2,file='DIFFEQ. OUT' ) 
open(unit=3,file='DIFFEQ. DAT' ) 


C READ INPUT PARAMETERS AND CONDUCT ERROR CHECKS. 
read(1,1000) numsys 


if((numsys. le. 0).or. (numsys. gt.4)) then 
write(*,1126) numsys 
stop 'The allowed values for numsys are 1 <= numsys <= 4. 
ena £ 
numplts = 2*numsys 
write(3,2000) numplts 


t 
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do 10 nprob=1, numsys 
read(1,1001) L, N, nstop, xsorce 


1£( CL. 190) ernie 128)) then 
write(*, 11245) probes Lene 
stop ‘The allowed values for “7G Sarre: §0)<= lel ee 
endif 


1f( CN: 1220) wormeaGh cc 128)) tee 
write(*,1124) nprob, a Nges 
stop 'The allowed values a ''N'' are: 0 <= N <= 128.' 
endif 


i£((nstop. lt. 0), ormtmstop. ciaole meme 
write(*,1127) nprob, nstop . 
stop 'The allowed values for nstop are: 0 <= nstop <= 850. 
endif 


if((xsorce. eq. ff Jeon xsoreencaqman se mchen 
Msorce.— fa 
elsei£((xsorce. equ Su or, (ssOrCeu camer mone 
xsorce = ‘S' 
else 
write(*,1128) nprob, xsorce 
stop ‘The allowed values for “ xsoree. sarc,mee or ueeoee 
endif 


C INITIALIZE EACH ARRAY TO ZERO BEFORE EACH RUN. 


dateawa/128*0, 0/34 1 29-0mG, 
data y/9/9"0. 07 yx 979-0307 


C READ THE COEFFICIENTS b(), a() AND THE INITIAL CONDITIONS 
C CyG Nee ie 


re2d(1,1002) (b(k), k=0,L) 
if(N. gt.0) then 
icstart = -N 
read(1,1002) (a(k)}, k=1,N) 
read( 1], 100ZISG@vCk)e k=i1cstrare.-1) 
endif 


C FOR xsorce = 'F' READ THE ARRAY x() FROM THE INPUT FILE. 
C FOR xsorce = 'S' CALL xgen TO GENERATE THE ARRAY x(). 


d £{(xsonee, cq. ft  yetnem 
read(1,1002) (x(k), k=0,nstop) 
else 
call xgen(x,nstop,nprob) 
endif 


C FOR TEST MODE ECHO INPUT PARAMETERS ONTO MONITOR (UNIT = *). 


if(mode. eq. 'Y') then 
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write(*,1007) 
weeee( =, 1120) nprob, infile 
weste(*,1110) 'INPUT' L', Ik 
write(*,1110) ‘OUTPUT’, 'N', N 
Pee *, 1112) nstcp 
write(*,1004) 
write(*,1005) (b(k) ,k=0,L) 
if(N.eq.0) then 
write(*,1131) 
else 
write(*,1006) 
write(*,1005) (Ca(Ck) ,k=1,N) 
endif 
write(*,1123) nprob 
pause 'END OF RUN, STRIKE <CR> WHEN READY TO CONTINUE. ' 
endif 


Peeve INPUT DATA INTO FILE: DIFFEQ. OUT. 


write(2,1008) ‘INPUT’, nprob 
wreite(2,1120) nprob, infile 
weate(2,1110) INPUT’, sit L 
piereec2,1110) ‘OUTPUT’, 'N', N 
veree (2,1112) nstop 
write(2,1004) 
pece( 2,1905) (bC(k) ,k=0,L) 
if(N.eq.0) then 
peeite( 2.1131) 
else 
write(2,1006) 
write(Z,1005) (CaCk) ,k=1,N) 
endif 


tee clll fHE INPUL SEQUENCE INTO FILE: DIFFEQ.DAT. 


peetee( 39,2001) nstop + 1 
write(3,*) ‘INPUT SEQUENCE x(ns) 
write(3,*) 'SAMPLE # (ns)' 
write(3,*) 'x(ns)' 

do 54 ns=0, nstop 


ii = ns 
werceco.2Ol0) 11, xCns) 
54 Continue 


CeCALE ditteq TO COMPUTE THE SQLUTION TO THE DIFFERENCE EQUATION. 
ea llvdutteaqt NG, 4a,b,x,yj;nstop) 
C WRITE RESULTS INTO FILE: DIFFEQ. DAT. 


WEreciorZOGl IN + nstop + 1 
write(3,*) 'OUTPUT SEQUENCE Y(ns)' 
Wieiter3,:) SAMPLE # (ns)’ 
write(3,*) 'y(ns)' 
doses eis— -N, NnStop 

ii = ns 
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write( 3,20 10) eae one 
55 continue 


C WRITE RESULTS INTO IE Sor EO: 


write(2,1008) ‘OUTPUT’, nprob 

WriteQ(2, biz) 

do 102 ns= -N, nstop 
wrate(2,1130)bns .x(is ya ycns) 

LOZ continue 

write(2,1123) nprob 


10 continue 


write(*,1121) . 
999 close(unit=1) 

close(unit=2) 

close(unit=3) 


if(ierm ct, Oy then 
Write: , 116) e1em: 
endif 


Cistevetesctees TNPUT FORMAT = vvtvedevedevest 


1000 format(il) 
LOOI Se formatC iene as el eee) 
1002 format(6f10. 0) 


PRET EN ATER GL OL SPR OD ES PL ON GL GD GL EU FU ETE OT OR ER OD EV IRR CS ENON OSES 


1004 format(t4, THE COEFFICIENTS b(O), bCL) 2. eC chee 

1005 a. formartCoC 1-676) 

1006 format(//,t4,' THE COEFFICIENTS a(1), ..., a(N) ARE:',/) 

L007 SformatG/7477777/// 77) 

1008 format(///,t16,a6,' DATA FOR PROBLEM # ',1i1,//) 

1110 format(t4,'THE NUMBER OF ',a6,' DELAYS: ',al,' = ',i3) 

1112, format(t4, THE VALUE OF nstep ilocos) 

1115 Oformat(1x,'DO YOU WISH TO RUN THIS PROGRAM IN TEST’, 
1' MODE ? CY) ND SG eee 

1116 Oformat(///,1x,'ERROR OPENING INPUT FILE, PROGRAM TERMINATED. ' 
1//, 1x ERROR CODE) Sian) aoe) 

1117 tormacGa) 

1118 Oformat(/////,1x,'TYPE THE NAME OF YOUR DATA FILE FOLLOWED’ , 
1' BY <CR>.',/,' IF YOU DESIRE TO MAKE A TEST RUN USING THE’, 
2' SAMPLE DATA ALREADY STORED',/,' IN THE FILE: DIFFEQ.TST , 
Sy TYPE: DIFFE@.TS@ <CR =) yeh ee) yee er) 

LISS torndecaiZ, 

1120 format(///,t4,'PROBLEM # ',i1,' | INPUT DATA SOURCEFILE: (‘,al2) 

1121 Oformat(//,' TABULAR OUTPUT DATA IS STORED IN FILE: DIFFEQ. OUT. ' 
1,/, | PLOTTING DATA IS STORED IN FILE: DIFFEQ. DAT. ') 

1122) Stormaciw.> 

1123 format(/yisgtc@ — )) Shimon PROBLEM {fee 16( - Dy /1) 

1124 Oformat(//, For problem #712, the valweutocs: ns eee 
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1'. This value is not allowed. ') 


DizoestOmmat(//, numsys = ',i4,'. This value is not allowed. ') 


17 . t 


1127 Oformat(//,' Fer problem #',i2,' the value for ''nstop 1S. ee 


173, . This value is not allowed.) 
ue 1? e t 


28 Oformat(//,' For problem #',i2,’ the value for ''xsorce'' is: ', 


lal,'. This value is not allowed.') 


Meo wetormat(té6é, ms ,ti6,'x(ns)',t35,'y(ns)') 

1130 format(t4,i4,t11,e14.6,t30,e14. 6) 

eeteetormat(/, THIS SYSTEM IS NON-REGURSIVE, I.E., N= 0.') 
2000 format(il) 

2001 format(i3) 

PUIO) format(el2.6,2x,e12. 6) 


end 


C SUBROUTINE: diffeq 
MecemrOskE: THIS SUBROUTINE COMPUTES THE SOLUTION TO A DIFFERENCE 
C EQUATIONS ALL PAKANBTERS DESCRIBING THE EQUATION AND 
C THE INPUT AND OUTPUT SEQUENCES x() AND y() ARE PASSED 
C TO THE SUBROUTINE BY THE MAIN®PROGRAM. 

Subroutine sd: feqiN,L,a,b,x,y,nstop) 

Boal Greece. cop), Vl -128instop),, a(l:N), bC0: L) 

do 500 ns=0Q, nstop 
y(ns) = 0.0 
foms0! k=0, max<:.,L) 

Pees =— (ns) aCk)~y(ns-k) + bCk)*x(ns-k) 

501 continue 
500 continue 

return 

end 
C SUBROUTINE: xgen 
PeeeereOsr: THIS SUBROUTINES ALLOWS THE USER TO GENERATE VALUES FOR 
C THE ARRAY x(). IF xsorce = 'S' THE MAIN PROGRAM WILL 
& CALL THIS SUBROUTINE. IF xsorce = ‘'F' THIS SUBROUTINE 
C WILL NOT BE CALLED BY THE MAIN PROGRAM. 


subroutine xgen(x,nstop,nprob) 
real x(-128: nstop) 


Sesscecss ont on os Vesescvesete veel lete les: woe eeasene eseseve ves: ale ws eveseveves's2 we tents ale ales Neve slesesloves'c s'e eves textes’ svevses'es’s testes vote sevedslctecs ae =": 
ze eur er fuew fy ee 


DEVELOP THE ALGORITHM FOR GENERATING VALUES OF oO my THiS SPACE 
THE STATEMENTS TYPED IN MUST FOLLOW STANDARD FORTRAN 77 RULES AND 
Revol FORTRAN 77 INTRINSIC FUNCTIONS SUCH AS: SIN(), COSC), .. 
NOTE THAT THE VALUE nprob CAN BE USED IN A LOGICAL 'IF' STATEMENT 
TO@MSTCH THE GENERATING FUNCTIONS TO THE CORRESPONDING SYSTEM 
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EQUATION READ FROM THE INPUT FILE IF MCRE Tian Ores eer 
EQUATIONS EXIST. AN EXAMPLE OF AN ALGORITHM GENERATING VALUES 
FOR x(a: 


isveve =O EXAMPZLE  worte 


if(nprob.eq.1) then 
do 1 k=0, nstop 
x(k) = 100.0 
il continue 
endif 
x(0) = 1.0 
do 1 k=l1,nstop 
x(k) = 0.0 
i continue 


MIMIAQqQIQIIgIaIaIgaqd 


Caededesedededededede dedegs feve We seve de te F886 Fe TENE Ge TE Te Te Te TE TET TE 8 TE TE ETE TE TELE TE TE TE TE TENE TE TET TETE TONE TEE ETE OTE TE TEN TES ERE IETS HE 


return 
end 


B. INPUT, OUTPUT DATA OF SECOND ORDER IIR NOTCH FILTER: P=0 


INPI> DATA FOR PROBLEM # 1 


PROBLEM # 1 INPUT DATA SOURCEFILE: diffeq.tst 
THE NUMBER OF INPUT DELAYS: L= 2 

THE NUMBER OF OUTPUT DELAYS: N= 2 

THE VALUE OF nstop IS: 400 

THE COEFFICIENTS b(0), b(1), ..., b(L) ARE: 


- SJOOQ0ETOO =. 3Z299S30R+00 200005200 


THE COEFFICIENTS a(1), ...) aCN}eARE: 
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more 2UETOO -. 980100E+00 


ieee 
how 


U 
WOMmMYWA NK WHF OF 


OUTPUT DATA FOR 


x(ns) 


- COOODDE+00 
. OOOCO0E+00 
- LOOOOCE+01 
- OOOOO00E+00 
- OOCOOODE+00 
-. OOOCCOE+00 
- OCOOOOCE+00 
. OOOC00E+00 
. OOOO00E+00 
. OOOOOCE+00 
- OOOCO0E+00 
- OOCOOCOE+00 
- OOOOCOE+00 
- COOOCOE+00 
- DOOCO0E+00 
. OCOCOCOCE+00 
. OCOOCOE+00 
. OOOCOCE+00 
. COCOOCE+00 
- COCOOCE+00 
. OOOOOOCE+00 
~eeoCece.C0 
- COOOONGE+00 
- OCOOO0E+00 
. DOOGOOE+00 
. COCOCCE+00 
. OOOOOCE+00 
2OOCUCOR 00 
. OOCOODE+00 
. OOOOOOE+00 
- OOOCOCE+O0 
COCCOCEs00 
. OOOOO00E+00 
-OOC0COE+00 
. COOCOOE+00 
- OOOCO0E+00 
- OOOCOCE+00 
. OOOOOCE+00 
. DOOOCCE+00 
~OCCCC 2100 
. COOOOCE+00 
. COOO0CE+00 
_OCC9005700 
. COOOO0E+00 
. COOOOCE+00 
. OOOOO0E+00 
. OCOOOOCE+00 
. OOOOCOE+00 


PROBUS 71 


y(ns) 


. OOOOO0E+00 
- GOOOOQ0E+00 
~IZOCCCE+OC 
sear) S603 10 
ei SOlZ 50.) 
-937446E-02 
lS eo 26 - Cn 
ele lE =O 
7 LONG 2E- 01 
oli Zeon =01 
waz £09 >- 02 
pele 2o 4b 0 1 
TIS OCE ~02 
By ores FESO 
77 Soe U2 
solic ole B= © 1 
eZ tee20h =O 
lec ZoL- Ui 
wile s068E-O01 
.494956E-02 
7 oOo 5 ee | 
Pes72 cn —03 
sleOegen-O 1] 
Wouea5 76-02 
leo £5 0n.—U | 
a ey SCr—() | 
SOM Se 0 1 
wiles Oo -0 | 
wos OUZOE —O0 2 
S02 30E=01 
eo voen -03 
. 158684E-01 
.452349E-02 
lov eulE-0 1 
767 66708-02 
ooo ele O02 
od Pose «O01 
fo coe ee R02 
2 owe cre 0) 
ee 2 So 02 
7 iS6/o@e-01 
O20 oe =O 2 
2S oar =0 1 
7 2 02 
966/655 -02 
s 1032265 -0 1 
. 606425E-02 
siz 136E-01 


i09 


. OOOOO0E+00 
. OOOODDE+00 
. OOOODDE+00 
. OOOOODE+00 
. OOOOO0E+00 
. OOOOOOE+00 
. OOOO0O0E+00 
. OOOO00E+00 
. OOOOD0E+00 
. OOOO00E+00 
. OOOOD0E+00 
. OOOODDE+00 
. OOOO00E+00 
. OOOODDE+00 
. OOOODDE+00 
. OOOOD0E+00 
. OOOOOOE+00 
. OOOOO00E+00 
. OOOOO0E+00 
. OOCOOOE+00 
. OOOOOOE+0G 
. OOOOD0E+00 
. OOCOD0E+00 
. OOOOO00E+00 
. OOO000E+00 
. OOOOO0E+00 
. OOOD00E+00 
. OOOOD0E+00 
. OOO000E+00 
. OOOOO0E+00 
. OOCOO0E+00 
. OOOO00E+00 
. OOO00CE+00 
. OOOOO0E+00 
. OOOOODE+00 
. O000V00ETOU 
. OOOOD0E+00 
. OOOOO0E+00 
. COCOO0E+00 
. OOOODDE+00 
-. OOOOCOE+00 
. OOOO00E+00 
. OOOOODE+00 
. OOOODOE+00 
. OOOO00E+00 
. QOOOOO0DE+00 
. OOOOD0E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. DOOODDE+00 
. OOOO00E+00 
. OOOOQ0E+00 
. OOOOO0E+00 
. OOOOOD0E+00 
. OOOOODE+00 


. 194032E-02 
1253248 =01 
; 225509 R Oe 
- IIS46ZE aoa 
. 599280288 =02 
9337 Sens O2 
. 695 9445<02 
, 619567 ba O2 
LOG so b=01 
p2a0072E -02 
.114349E-01 
pele = 02 
Oey te 0 1 
-484919E-02 
solo USE -O02 
wi POPS On=02 
N625577 2202 
pole sk 02 
. 294043E-02 
elOe 07 Ol 
soe OoUa=O02 
OOO S 7 E=01 
To US SOE OZ 
mop o7 OSE. O2 
SOD Ome OZ 
pote aer OZ 
$eaot son -02 
ol 7 a2 9b -02 
. 940264E-02 
| Gee st 8) 6: 
Poe UeOE=U2 
NETS) OWA 
.809430E-02 
Ue pe fs So 0/4 
mOUO DO? ES0Z 
. 746086E-02 
725d one = 02 
.847247E-02 
NOs eson 03 
nop Leou E02 
e277 Z2O5E-02 
ii eo226R-02 
-464753E-02 
»-594701E-02 
.651714E-02 
5 Sheet hd S10) 
7 OU LOBE -—O2 
. 109744E-02 
ay Ola 02 
Sleoly ShsvZ 
4 J MSOEK es (042 
Nels LS) We S10 M2 
pov 24-02 
. 565640E-02 
mofo eg -OZ 
7O7 OGD 5h O2 
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135 
136 
137 


138 


. QOOO00E+00 
- QOOOO00E+00 
. DOOOQ0E+00 
. OOOO00E+00 
. DOOO00E+00 
. QOOQ00E+00 
. QOOOQOO0E+00 
. OOOOO00E+00 
. QOOQOQ0E+00 
. QOOO00E+00 
. DOOO00E+00 
. QOOO00E+00 
. QOOOO0E+00 
. QOOO00E+00 
- QOOO00E+00 
. QOQOQQ0E+00 
. QOOOO0E+00 
. QDOOO00E+00 
- OOO900E+00 
. QOOGCGE+60 
. 000000r+90 
- 000000£4+00 
. OOOD009E+00 
. QOOOO00E+00 
. QOOOQDE+00 
. QOOO00E+00 
. O0O0000E+00 
. DOOOQ0GE+00 
. 0000002£+00 
. QOOO00E+00 
. QOOQQODE+00 
. OOOODDE+00 
. OOOOO00E=00 
. O0O00C0E+00 
. DOOOO0E+00 
. QOOO00E+00 
. QOOOCOE+90 
. JOOOD0E+00 
. OCOOCOE+00 
. O0O0090E+00 
. QOOOOCE+00 
. OOOCOCE+GO 
. QOOOQ0E+00 
. QOOOO00E+00 
. OOOOCOE+00 
- OOOOO0E+00 
- 0O00000L+00 
- DOOOQ0E+00 
-. QOOO00E+00 
. OOOOODE+00 
. QOOOO0E+00 
. QOOOO00E+00 
. QOO000E+00 
. QOOOQ0E+00 
. QOOD00E+00 
. QOOOQ0E+00 


, Hes 7G1E-02 
.713443E-02 
. 925062 -03 
5 6057 252-02 
ole gore O02 
. S27 ee 18 2 
. B74 SF =) 2 
- SB08S5E-02 
. 603424E-02 
. 17421 Ma-02 
- 648895E-02 
. 423535958 -03 
. O21 663E-02 
2 SOO O02 
» 9276245 -02 
.416741E-02 
. 37 eZ IE O02 
wee 9 E-O2 
ge oe =02 
poe eee -02 
. 182246E-04 
.5754442-02 
Io lG22n-O02 
.500764E-02 
. 3530456E-02 
nog aol -O2 
wees on =02 
wet ious -02 
ToSOL0SE-OZ 
~o20645n-03 
s DOU SIGE -O2 
. 16260702 
24727262 -02 
eco oe -O2 
7 20500e8-02 
. 410886E-02 
2222906b=02 
. 4760592-02 
. 6140588-03 
. 486845E-02 
. 100441E-02 
.444019E-02 
. 244937E-02 
. 3543708 -02 
; S567 ems -02 
. 229540E-02 
4250098 -02 
.845509E-03 
sao Og55-02 
Ce 76225-03 
. 8 50762-02 
. S03 5E -02 
. se0950e-02 
suet Ze 0 2 
POAC LL 15) gl OV 
movo72/h-O2 


11} 
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. OOOO000E+00 
. OOOO00E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. O00 CRE 00 
-. QOOOO00E+00 
. OOOO000E+00 
5 GLONONONOE Do of O10, 
. QOOOOD0E+00 
. OOOOO00E+00 
. QOOOO0E+00 
. QOOOOOD0E+00 
. QOOOOQ00E+00 
. QOOOOQO0E+00 
- OOCOOC0E+T00 
- DOOOCCETOO 
. OOOOD0E+00 
. OOOO00E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. QOOO00E+00 
. 000000800 
. OOOO00E+00 
. QOOOO00E+00 
. DOOO00E+00 
. QOOOO0E+00 
. OOOO00E+00 
. DOOOO0DE+00 
. OOOOO0DE+00 
. DOOOO0E+00 
. QOOOO00E+00 
. QOOOO00E+00 
. OOOOO0E+00 
- OOOO00E+00 
- OOOOD00E+00 
- OOQ000ETO0 
. OOOO00E+00 
000000800 
. OOOOODE+00 
- OOOOQ0E+00 
. OOOO00E+00 
. OOOO00E+00 
- OOOOODE+00 
. OOOO00E+00 
- DOOO000E+00 
. OOOOO0E+00 
- OOOOD0E+00 
- QOOOO00E+00 
 COCO00ERD 0 
. OOOOO0E+00 
- OOOO00E+00 
. QOOOO0E+00 
- QOOO00E+00 
0000005 --00 
- OOOO000E+00 
. OOOOO0E+00 


. L029 ter —02 
-405146E-02 
JoeZGUtor-03 
» S507 elE O02 
lS yore OZ 
Bes Pees) e) 2 RO 
. 263944E-02 
. 23235-0902 
» 3353 5m = 02 
2 117087E=02 
7507 206E=02 
.642924E-04 
. So evie O02 
e227 en - 02 
9097 2 1ESO2 
. 224087E-02 
227024r,-02 
oe) 5)518)/ 10 04 
, LZ Ss See 02 
| OSS eS Oey 
wsoonlOn—03 
re Oe 2 
9478-03 
27 270, -02 
soos TUR =O02 
eeeo7 Can-O2 
Poor 0b -02 
plea I50H=02 
s2eo0glE 0 2 
Poo JUS 
. 303344E-02 
POOG> 4-03 
ae oe OR 02 
5 LS ey 1S tess O rs 
weeolg2n=O2 
peo or, O2 
. 139544E-02 
~oeon or 02 
rae? o 10-03 
wo oa 9r.-02 
Too Jot — 03 
e2o vo or, O02 
mie IOZE=02 
5a UO 5) Se O re 
5 MSE SS 
wee or -02 
ee Oar 02 
FOU 59 a7 - 0 
72 SASH OSC 8 
> 2405128-03 
Peodgoy O02 
e202 
e 201267 8-02 
ploy Vest 02 
leo ee =O 2 
molLOGbSE 02 


214 
215 
216 
Za 
218 
fasts, 
220 
Zou 
222 
223 
224 
225 
226 
227 
228 
229 
230 
ZS 1 
Zoe 
253 
234 
oo 
230 
or 
220 
ISS, 
240 
241 
242 
243 
244 


245 


264 


. OOOOO0E+00 
. OOOOODE+00 
. QOOOO00E+00 
. OOOO00E+00 
. QDOOOO0E+00 
. OOOO00E+00 
- OOO00CE+00 
. 00000C24+00 
. OOOO00E+00 
. OOOO00E+00 
OGROCCErOO 
. OOOOO0E+00 
. OOOOO00E+00 
. OOOOO0E+00 
. OOOOO00E+00 
. OOOOO0E+00 
. QOOO00E+00 
. OOOOO0E+00 
. OOOO00CE+00 
. OOOOO00E+00 
. QCOOOO0E+00 
. QOOOO0E+00 
. OOOOO0E+00 
. OOCOO0E+00 
. OOOOOCE+00 
. QOOOO00E+00 
. OOOOO0E+00 
. OOOO00E+00 
. OOO000E+00 
. OOOO00E+00 
. OOOO00E+00 
. OOCOO0E+00 
. OOOCOO0E+00 
. OOOOD0E+00 
. OOOO00E+00 
. OOOOO0E+00 
. OOCOOOE+00 
. OOOO00E+00 
. QOO0COE+00 
. QOOO00E+00 
. OOOOOOE+00 
. OOOO00E+00 
. OCOOO00E+00 
. OCOOO0E+00 
. OOOO00E+00 
. OOOO00E+00 
. OOOOO0E+00 
. OOOO00E+00 
. OCOOO00E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOD0E+00 
. OOOO00E+00 
. OOOOO0E+00 
. OOOOQDE+00 
. OOOOO0E+00 


PO ck 05 
22 UGE = 02 
7/308 255-0% 
PAPEL ARADO 
ou. /IoE-03 
> eS Can - 02 
- 142160E-02 
- 140794E-02 
eer eob-02 
; 7eeIooE-03 
207 Sot 2 
- 674610E-04 
5 AUS YES 0) R10) 74 
Toutes 5b. -03 
pede = 0 2 
ele sok - 02 
mis oOIOR<-02 
wlG 23255 -02 
Jota oE=03 
wiGoGizE-O2 
SO SOIE = 03 
. 188968E-02 
.443136E-03 
: O58 7E-02 
> ge 137-03 
Tle ar 02 
lo 20 k= 02 
- 847595E-03 
Gale Si-O 2 
pe o oT E03 
Tl S0DSE-O02 
we SOen, Uo 
. 159804E-02 
poueoe7 E05 
ele iober 02 
2 W2e003E-02 
.864516E-03 
1490755 -02 
Ho] 499E-03 
ey oooh OZ 
yee ona 0S 
| L490 10E-—02 
.660574E-03 
eee olLE-O02 
Oo 02 
eee Z2ob-03 
oe oO = 02 
.415134E-03 
. 143365E-02 
. 661320E-04 
leo 00E-02 
ee 20s 0 5 
. SS orE=02 
glee =O 3 
HeoweozZE-03 
wilson 07 


i 


270 
27a 
Zi 
275 
274 
Zi 
Pidils 
pag | 
278 
279 
280 
7) 35 
ZE2 
269 
284 
235 
286 
ZO7 
288 
280 
Z70 
Zu 
(ahs i 
yi 
294 
Paes |, 
296 
297 


. OCQO00E+00 
. DOOOOOE+00 
. OOOO00E+00 
. OOOOO0E+00 
- VOOCOCE+00 
. DOOO000E+C0 
. OOOOO00E+00 
. OOOOO0E+00 
. QOOOOO0E+00 
. DOOOO0E+00 
. QOOOOD0E+00 
. QOOOO00E+00 
. OOOOO0E+00 
. DOOO00E+00 
. OOOOO0E+00 
. QOOO00E+00 
. QOOOO0E+00 
. O0O000'E+00 
. 0OCO0CETC0 
. OCOOOO0E+00 
. QOOOOO0E+00 
. COCR COer00 
. QOOOO0E+00 
. OOOOO0E+00 
. QOOOO0E+00 
. DOOOO00E+00 
. QOOOOO0E+00 
. QOOOOO0E+00 
. OOOOO0E+00 
. QOOOO0E+00 
. QOOOOO0E+00 
. QOOOOQ0E+00 
. QCOOOOD0E+00 
. DOOOO00E+00 
. DOOOCOE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOO0E+00 
- QOOOOD0E+00 
. OCOOOO00E+00 
. DOOOO0E+00 
. QOOOOOOE+00 
. QOOOOO0E+00 
. OOOOOOE+00 
. OOOOO0E+00 
. QOOOO0E+00 
. OOOOO0E+00 
. OCOOOO0E+00 
. OOOOO0E+00 
. OCOOOO0E+00 
. OOOOODE+00 
. OOOOO00E+00 
. OOOOOD0E+00 
. QOOOO0E+00 
. DOOOODE+00 
- OO0CCOEL00 


. 4000 7ZE=035 
, 1296775 -0 2 
. 202 (OTE 04 
, I2766er- 02 
. 399054E-03 
LIZ Ser 02 
-/ 611408 -03 
. 848092E-03 
. LOZSS0E 02 
49203 
, ero 7-02 
.976162E-04 
lilo o)3 04 
29263 2n=035 
. LOds697 B=O02 
,6355545-03 
. 8202 30h =03 
. 895491E-03 
- SIG sa68-03 
, LOgye7E-02 
rie ss05b-03 
On Ser = 02 
2005 505.-03 
oo TZ E03} 
ico Sow 0 5 
. 798744E-03 
7) 7O6o7 E05 
2 SGe 7 E05 
po oao 7 LE = O03 
. 207643E-03 
. 984873E-03 
wien 03 
= 7 E03 
.424267E-03 
HiOee2Z1E-03 
. 668821E-03 
. S3083 78-035 
83007 0E "05 
e240 207h-—05 
peste hay os)! {Lah ag O Re, 
. 540483E-04 
| 6297 203 
~9566275-03 
J7 ole OR =035 
7 A= OS 
192052 5E=035 
. 734302E-03 
a2 5) Se S 
. 810663E-03 
127 5625 E-05 
Be s5 109,105: 
3257 7 OE =03 
peo e S15) al ON, 
-483492E-03 
a 20CORE=05 
- 645634E-03 
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326 
Sy 
626 
O29 
550 
eo! 
Dae 
oa. 
334 
255 
220 
JO7 
338 
oo 
340 
341 
342 
343 
344 
345 
346 
347 
348 
34S 
520) 
seul 
Bo 
Pome) 
354 
S155: 
o50 
307 
326 
eo, 
360 
Sie! 
SoZ 
363 
364 
365 
366 
367 
368 
369 
270 
S) pl 


573 
374 
SYS, 
378 
5) 717 
378 
379 
380 
Sou! 


POUCOCCESOC 
. OOOOQODE+00 
. QOOOO0E+00 
. OOOOO0E+00 
-. OOOOO0E+00 
. OOOOO0E+00 
. QOOOO0E+00 
. QOOOO00E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOO00E+00 
. OOOO00E+00 
. OOOOO0DE+00 
. OOOOQO0E+00 
. QOOO00E+00 
. OOOOQ0E+00 
. OOOOO0E+00 
. OOOOOOE+00 
. QOOOO0E+00 
. OOOOOCOE+00 
. OOOOQ00E+00 
. QOOOD0E+00 
. QOOOOOE+00 
. OOOOO0E+00 
. DOOOQO0E+00 
. OOODOCE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOO00E+00 
. 00000Cc+00 
. OOOQOO00E+00 
. OOOOO0E+00 
SoCs Orn) 
. OOOOO00E+00 
. OOOOOOE+00 
. QOOOOOE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. QOOOOO0E+00 
. QOOOO0E+00 
. QCOOOO0E+00 
. QOOOOO0E+00 
. OOOOO00E+00 
. OOO000E+00 
. QOOOO0E+00 
. DOOOOQ0E+00 
. OOOOO00E+00 
. DOOOO0E+00 
. COOOODE+00 
. GOCCOOE+T00 
. OCOOOOOE+00 
. DOOODDE+00 
- COUOCOETOO 
. OOOO00E+00 
. QOOOO0E+00 
. OOOOO0E+00 


$2992 535-05 
27 Noe 2-02 
UTS 765-04 
So 2007-05 
lo ZO CDE 0S 
- 654736E-03 
-404790E-03 
so 0G Ion -O3 
. 564390E-03 
7 Sloss =05 
. 656042E-03 
.891797E-04 
won ouLOk=O 3 
. 134443E-03 
weteoy 2E-03 
5 B43 Oe US: 
.492056E-03 
PaS0255b-03 
oe on O03 
- 586244E-03 
S27 Ter 03 
.614406E-03 
. 843958E-04 
. 974335E-03 
A PAU TAD AOE: 
.473097E-03 
Ped eae): 
. 324162E-03 
po2IG1GE-05 
. 145681E-03 
pooeTO2E-03 
.416829E-04 
,oee226r-03 
ave OR 09 
7 om sole 05 
Poets 0SE-03 
. 323564E-03 
»461457E-03 
. 164877E-03 
5 EOE 7! ail: 
sD57 0535505 
-494751E-03 
ISHS SIS SIO): 
.429248E-03 
, S00958E—03 
7 ogee 2h = 02 
~tU0CZ59R-03 
Be ac OMe) ocd ORG, 
poe) 03 
. 246287E-04 
.456254E-03 
J1Z639SR 2035 
-405474E-03 
. 257656n-03 
Ze OL-0s 
72D 5 I6E 05 
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oo2 . QOOOO0E+00 -IS8Go7Rh=05 
Biss . OCOOOO0E+00 -410831E-03 
384 . OOOOO00E+00 P4953 56704 
385 . OOOO000E+00 .419005E-03 
386 . OOOOQ0E+00 .896748E-04 
387 . QOOOQ0E+00 eno Sul=03 
388 . OOOOO0E+00 2 oz =03 
309 . QOOOO00E+00 , SUSGI E03 
390 . QOOOO00E+00 2309 5¢5R—03 
oi . QOOOO00E+00 21947 25b—03 
SoZ . OOOOO0E+00 ~367499E-03 
Chr). . QOOOQ0E+00 .697972E-04 
394 . OOOOO0E+00 sood2 Fak -03 
595 . QOOOOOOE+00 . 580256E-04 
BO . OOOOO0E+00 .356444E-03 
Boy . OOOO00E+00 .174472E-03 
58 . OOOOODE+00 a2 ee O 3 
oo? . OOOO00E+00 n2e 27 OR) 3 
400 . OOOOD0E+00 le SOCh=03 


woo een ane === == END OF PROBLEM ## 1 ---------------- 


INPUT, OUTPUT DATA OF THE AUGMENTED IIR NOTCH FILTER: P= 1 


INPUT DATA FOR PRGSUEHS 1 


PROBLEM # 1  IN:UT DATA SOURCEFILE: diffeq.tst 
THE NUMBER OF INPUT DELAYS: L= 3 

THE NUMBER OF OUTPUT DELAYS: N= 3 

THE VALUE OF nstop IS: 400 


THE COEFFICIENTS b(0), b(1), ., bCL) ARE: 


. SJD000E+T00 =.329900E-02 .8SIeernOGee, 326631700 


THE COEFFICIENTS a(1), ..., a(N) ARE: 


. OOOOOOE+00 -.871246E+00 -. 323364E+00 


OUTPUT DATA FOR PROBLEM # 1 
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x(ns) 


. QOOO00E+00 
. OOOOO0E+00 
- QOOOODE+00 
, LOGUE OETO 1 
. OOOOO00E+00 
. OOOOQ00E+00 
. OOOOQ0E+00 
. OOOOO00E+00 
. DOOO00E+00 
. OOOO00E+00 
. OOOO00E+00 
. QOOOD0E+00 
- OOOOO0E+00 
. OOOO00E+00 
. OOOO00E+00 
. QCOOOO0E+00 
. OOOOQ00E+00 
. OOOOODE+00 
. OOOOD0E+00 
. DOOOQ0E+00 
. OOOOODE+00 
. OOOOQ0E+00 
. OOOO0CE+00 
. OOOODDE+00 
. OOOO00E+00 
. OOCOOOE+00 
. OOCCOOE+00 
. QOOOO00E+00 
. OOOOQ00E+00 
. OOOO00CE+00 
. QOOO00E+00 
. OOOOOOE+00 
. OOOCO0E+00 
. OOOOO0E+00 
. OOOCOOE+00 
. OOOO00E+00 
. QOOO00E+00 
. OOOOQ0E+00 
. QOOOOOOE+00 
. OOOO00E+00 
. DOOO00E+00 
. OOOOO00E+00 
. OOOOO00E+00 
. QOOO00E+00 
. OOOON0E+00 
. QCOOOD0E+00 
. OOO000E+00 
. OOOO0OCE+00 
. O0000C=E+00 
. DOO000E+00 
. QOOOO0E+00 
. 000080L+00 
. OOOOD0E+00 
. QOOO00E+00 


y(ns) 


. OOOO00E+00 
. OOOO000E+00 
. OOOOODE+00 
- JABV00E+00 
Poe OUn =O 2 
1 Obsead E = O4 
aS 74C0r 02 
15 7469cE-01 
- 141862E-01 
wiGme 75-01 
ei 25058501 
2 Yale 5 — 0 2 
Plos7z4aZ2b-01 
. 185964E-02 
eo oo O0r —0 1 
ee s598-=02 
7ia5 L13E-01 
mio We) | 
SO ZOZ9 <0 | 
) US SUIS IB eM 
.494966E-02 
Ooo so8 =) | 
meow oor =-035 
. 160940E-01 
novel E02 
loo 56r <0 1 
06 7475-01 
pe Onwey 9 b= 01 
5 SEO eal 
5 EOS SNe 
lau s92b=0 1 
ooo” 6-03 
wie t- 0 | 
» 452347E-02 
2ioOS02E—01 
more g00r—O2 
Io Sos Eh - U2 
pilesZ 2b Ot 
Boo 25 258-02 
Poo oo SE -01 
sIZ2e268E-02 
S230 75GE=01 
moog Sele O02 
= NAS I SISO 
./29410E-02 
BICo7 77 2-02 
mle Ss son -0i1 
. 606436E-02 
le Sos 01. 
- L94O41E-02 
7 £25324E-01 
~22s500ZE=02 
iiS4o5F-01 
Pog 7cU5E-O2 


Dele 


. OOOOO0E+00 
 OOCUGOES00 
. QOOOOQO0E+00 
. DOOOO0E+00 
. QOOOO0E+00 
. QOOOD00E+00 
. OOOOODE+00 
. QOOOOOOE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OCOOO00E+00 
. DOOOO0D0E+00 
. QOOOOOE+00 
. OOOOO0E+00 
. OOOO00E+00 
. OOOO00E+00 
. OOOOO0E+00 
. DOOO0E+00 
. QOOOOO00E+00 
. OOOOODE+00 
. QOOOOO0E+00 
OCUOOCE 00 
-. QCOOODOD0E+00 
- OOOOOOE+00 
. OOOOO0E+00 
. DOD00E+00 
- OOOOQ0E+00 
. DOOOOO0E+00 
. OOOOODE+00 
. OOOOO0E+00 
. OOOO00E+00 
. QOOOOGOE+00 
. OOOOO0E+00 
- OOOOQ0E+00 
. QOOOOO0E+00 
. DOO000E+00 
. DOOOOGE+00 
. DOOO90E+00 
. OOOCOOE+00 
. QOOOOO0E+00 
- OOOOO0E+00 
. DOOOOQO0E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOO00E+00 
- QOOOOO0E+00 
. QCOOOODE+90 
. OOOO00E+00 
. OOOOO0E+00 
. QOOOO0E+00 
. OOOOO0E+00 
- OOOOO00E+00 
. QOOOOODE+00 
- QOOOO0E+00 
- QOOOOO0E+00 


. $og7S50E-02 
. 895943E-02 
,619s79-02 
, 10e23er-01 
~ 2oU0@en-02 
. LT4S50R<01 
. 1320esE =02 
 LOGV a =O 1 
-484913E-02 
wOgo? TOE a2 
TORS S11 
neO2Zsco on -O2 
Jo ltevE-o2 
. 294052E-02 


- J63500E—02 
sOoo FEBS OZ 
7Oo7 Goon 02 
-olgg22 8-02 
S49 1508-02 
,o27oock-O2 
Poa Oem r= (2 
O07 0OVE =03 
29 250e2E=02 
ee ean? 
»S094S46E-02 
Zoom —OZ 
POUyM Ion -O02 
macus TE-02 
ae cre = OZ 
.847244E-02 
nen Og oi -03 
65 LOO0R-02 
5 A EN eNO 
Os ZAR 02 
-464744E-02 
oo O02 
Ooty a7 E=02 
FOC 7 So4E-02 
. 760104 -02 
SOE 7 SOs 
7] See = 02 
ls O leo OZ 
ICU Ee =02 
- 383442E-02 
57 Soe, 02 
TOs zr 02 
mov Zoo OZ 
. 678848E-02 
-145789E-02 
o7 13539502 
POZE 767 E03 
70057 24E-02 
yo zoek O02 
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- QCOOOO0E+00 
- QOOOOCE+00 
. OOOOO0E+00 
- OOOOO0E+00 
. QOOOO0E+00 
- OCOOOO0E+00 
sOCCOUCk.00 
- OOOOO0E+00 
. OOOOOD0E+00 
. OCOOOOOE+00 
. OOOOOOE+00 
. OCOOO00E+00 
. OOOOO0E+00 
. OOOOOOE+00 
. COOOO0E+00 
. OOOOOOE+00 
. COOOO0E+00 
. OOOOO0E+00 
. OOOOOOE+CO 
. OCOOOOOE+00 
. QCOOOO0E+00 
. OOOOOOE+00 
. OOOOODE+00 
OCCCCOE.00 
. OCOOOODE+00 
. OCOOOOCE+00 
. OOOOOOE+00 
. OOOOO0E+00 
/OCCOCOETO0 
. OCOOOOOE+00 
- COOOOCE+00 
. OOOOOOE+00 
. COOOODE+00 
. OCOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. DOOOOOE+00 
. OCOO00E+00 
. OOOOOOE+00 
. COOOOOE+00 
. OCOOOOOE+00 
. OCOOOOOE+00 
. OOOO00E+00 
. COOOOOE+00 
. OOOOO0E+00 
. OOOO00E+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOO0E+00 
. DOOOOD0E+00 
. OOOOGOE+00 
. OOOOOOE+00 
. OCOOOOOE+00 


. S27 ok =02 
-487450E-02 
- SOU SIGE <02 
. 603416E-02 
eee 2 5-02 
. 648890E-02 
o9eeOecr-O3 
Oeenconn = 02 
gE 7) Sj) 31 B16) 7 
ao 2) Oi 02 
4167 318-02 
Jo” eee =O2 
Paes 7 LE -O2 
~ 19S993E-O2 
i2o7 7 On-U2 
. 181364E-04 
.575441E-02 
Ig egor =02 
5 SC Sh 10 
wes 2080-02 
S91 S22E,-02 
.469510E-02 
oa 67 E =02 
5 2 SO Ue eIG ys 
EVA PL eNOS 
2 OSU 07 
le 2 9 5-02 
~472725E-02 
Pear O2 
Bylo, 2 ch sio) dale V4 
een 2B = OZ 
mee eo lor -O02 
Ta OUSZE (2 
Poles sh=03 
-486840E-02 
, LOCSSZE=02 
-444016E-02 
~2UG972 7 GeO 2 
2259457 UE S02 
Spe Se) tig wi 0 
. 229544E-02 
-425600E-02 
504557 78-03 
waa 502 9h =02 
| Segoe) 13 20s 
-415072E-02 
SIG or = b= 02 
. 340950E-02 
soUBTa2E-02 
5 PE PAS 1013) 3, 10% 
1 eer IE = 02 
O22 2h S02 
3405 I39h-02 
5 SVE S) S) 2) dint Oi 
Be eey45) (Noes 
re 2 2.) OP 
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fee 
164 
165 
166 
ie 
168 
ie? 
170 
a 
gy 
7S 
174 
175 
176 
lay 
1 
179 
180 
teal 
132 
183 
184 
185 
186 
187 
188 
189 
190 
Ieeyal 
122 
193 
194 
eS: 
196 
MS) 


ie 
200 
200 
202 
PAGS) 
204 
205 
206 
207 


209 
Zn MAG, 
ZN Al 
A? 
ZL) 
214 
ANS 
ZL 
AMG 
218 


. OOCCCOE+00 
. OOOOQ0E+00 
- QO0NNCL+00 
- OOOOO00E+00 
. OOOOOOE+00 
. OOOO00E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. QOOOO0E+00 
. OOOOO00E+00 
. OOOQO00E+00 
. QOOO00E+00 
. QOOODDE+00 
. OOOOO0E+00 
. QOCOCOETO0 
. OOOOO0E+00 
. OOOO00E+00 
. OOOOO0E+00 
. QOOO00E+00 
. OOOO00E+00 
. OOOO00E+00 
. QOOOO0E+00 
. OOOO00E+00 
- QOOO00E+00 
. OCQOOO0E+00 
. QOOO00E+00 
. OOOO00E+00 
- OOOO00E+00 
- OOOOO0E+00 
. QOOOO0E+00 
. OOOOO00E+00 
-. OOOQ00E+00 
. QOOOO00E+00 
. OQOOOO0E+00 
. OOOOD0E+00 
. JO0O0OE=00 
. QOOOD0E+00 
- QOO090E+00 
. OOGO00E+00 
. 000000£+00 
. OOOOQ0E+00 
. OOOO00E+00 
-. QOOO00E+00 
. OOOOODE+00 
-. QOOOO0E+00 
- QOOOOO00E+00 
- QOOOO0E+00 
. OOOOO0E+00 
. OOOO00E+00 
. OOOOO00E+00 
- OOOOODE+00 
. QOOOO0E+00 
. QOOOO0E+00 
- OCOOO00E+00 
. QOOOO0CE+00 
- OOOO00E+00 


200 2A O 2 
20599 an= 0 2 
~2ocoooe 02 
,e35e42E-02 
Ue Zz = 02 
, S07 SOE = 02 
. 642165E-04 
apo) fe Ti oN AO 
Be AES 5) 7/10, AOA 
7097 Tome OZ 
»224077E-02 
S22 702050 2 
. 22557 0n-02 
tee Ze 02 
, SOLOUOE 0 2 
5657-03 
o> ONG a0 2 
Be SS (eS lod 10S: 
o 2927208 =02 
icp OUr=02 
2267 698-02 
502 Or er =O 2 
els yes 1230, 103 
~298063E-02 
e077 OE-03 
Cos eG 02 
. 668474E-03 
nee ceon-O2 
Ga 2 i -02 
acer O2 
2252 0OE 02 
logo e7 E02 
mcoe) G5 E-02 
»487564E-03 
cei oeee=O2 
»437834E-03 
- 2o7 5 7en-02 
oes = OZ 
weep! O02 
Peigey de EsO2 
- 1a1S25E=02 
-23/63565-02 
. OUS0USE-03 
MAS) 72/85 5) 210.100 
- 240445E-03 
SIDES Se 0s 
OZ 727 = 02 
720 6CE-—02 
s67 Uso = 02 
» 142134E-02 
A PODS S10 
sOUs0s 7-05 
eee = OZ 
i202 IS =04 
oe ey OZ 
vewoy er. =0 > 
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274 


. OOOO00E+00 
. OOOOO0E+00 
. OOOOO0E+00 
- OOOOO0E+00 
. OOOOOCE+00 
- OOOOOOE+00 
- OOOOOCE+00 
. OOOCO0E+00 
. OOOOO0E+00 
. OO000CE+00 
. 0O0000Cz+00 
. OOOOO0E+00 
. OOOOO0E+00 
POV Cr O0 
. OOOOO0E+00 
. OOOOO0E+00 
. QOOOOO0E+00 
. OCOOOO0E+00 
. OOOOO0E+00 
- OOOOCOE+00 
. OOOOO0E+00 
. OOOOOCE+00 
. OOOO00E+00 
- OOOO0DE+00 
- DOOQO0E+00 
. OCOOO00E+00 
. OOOOODE+00 
. OOOCO0E+00 
. OOOOCOE+00 
. OOOOOO0E+00 
- OOOOODCE+O0 
. OOOOO0E+00 
. OOOOO0E+00 
- OCOOOO0E+00 
. DOOOOCE+00 
. OOOOO0E+00 
. OOOGOOE+00 
. OOOCOOE+00 
. OOO00CE+09 
. OOOOOCE+S0 
. OOOO00E+00 
. OOOOO0E+00 
. OOOOOOE+00 
. OOOOO00E+00 
. OOOOO0E+00 
. OOOOODE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOCOOE+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOOODE+00 
. OOOOCOE+00 
. OOOOD0E+00 
. OOOCOOE+00 


Pig SOChs02 
G2 e202 
le O79 SE 02 
ear u- 02 
e700 =03 
sou 7 Soar -O02 
pemeta2 OG 
. 205485E-02 
Ode Sob -03 
eZ te - O02 
PIs 7esk<02 
wiooOlon 02 
ioe 7 E02 
POM ooE-=03 
. 186605E-02 
. 184036E-03 
mlSo Osh 02 
.443071E-03 
.170584E-02 
~ 227 000E=03 
.134294E-02 
lee OS OE -—02 
76470138 -03 
ey ose =02 
. 9 207 E-03 
ao le05ar-O2 
weo7 [0G -03 
tome Oth. = 02 
- §18618E-03 
el PaE) 6 NP) F182 
Piz Sy Or 0 2 
. 864528E-03 
meg OLE O2 
Ho 9495.,-03 
» LoWSS3E-02 
si Z521E-03 
Pe 7007E-02 
"eoUs0E-03 
elZa2506-02 
71057 308-02 
Psiloc Soo Os: 
wee ego - 02 
eevee: =O 5 
. 143360E-02 
. 660854E-04 
sles geok-02 
“De eeE-05 
peligeye E-O2 
7 oOdsa7 03 
~OO255cGr,-03 
Pitesogr 02 
~200295n-03 
mlzge7 1-02 
25 1AZE-04 
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D. INPUT, OUTPUT DATA OF THE AUGMENTED HR NOTCH FILTER: P=2 


INPUT DATA FOR PROBLEM # 1 


PROBLEM # 1 INPUT DATA SOURCEFILE: diffeq.tst 
THE NUMBER OF INPUT DELAYS: L= 4 
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. OOOOQ0E+00 
. OOOOO0E+00 
. OO0000E+00 
. OOOOO0E+00 
. OOOOO0E+00 
. OOOO00E+00 
-. OOOO00E+00 


woo 74E,4+09 
IOC lOUEO9 
TOo222509 
2 2205 OEeO9 
00 79 0be09 
696671709 
- LOGCSOET10 
eS or 10 
zoe 7 ee O 
Pig 765% LO 
. 153410E+10 
a O07 20be 10 
Se oye kej 5 omaae 
. Ze ZOE 10 
. 23527 6b 0 
p2OUsz abs) O 
G28) Sieloleae 0, 
. 324243E+10 
, S608 29R10 
-401544E+10 
-446852E+10 
16972 75b LO 
. 553384E+10 
. 6156255710 
/ 6658078710 
. 762640E+10 
. 848693E+10 
. 944456E+10 
Ose Zz Ed 
Pre 262 bell 
eels 0) 1592S IE 
. 144846E+11 
. He PORTA 
1937 CE 
SLE ie Walang 
222 eee tl 
.247208E+11 
‘27 oe ier 1 
-S0U61S3ET11 
. 340687E+11 
Beli. cildan Ika 
.421908E+11 
.469514E+11 
Sa2Z2692E eA | 
-581448E+11 
647050711 
~A20007 Ea 
oie C7E I. 
§Og leet 11 
wJoZooegk+ 11 
ay IMI O45 5612 alleys 
wie ZeooE 12 
Soy OO 1 2 
woe Ze 12 
7 169364E+12 
. 1O84a5E+12 


1S! 


eo . OOOOOOE+00 . 209741E4+12 
o9:2 . COOOOOE+00 -. 233408E+12 
SS): - COUCOUE+00 He Geel 2 
394 . QOOOO0E+00 ee AS SSS) ia aes 
BIS . OOOOOCE+00 so 2oO9r +12 
396 . COOOOOE+00 =. Oo aoa | 2 
O97 . QOOOOO00E+00 -3982s0ud 1 2 
398 . OCOOOOOE+00 -. 443305E+12 
SEN, . QOOOQO0E+00 - 493325E+12 
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wo nnn nee --- eee END OF PROBLEM # 1 ---------------- 


INPUT, OUTPUT DATA OF THE AUGMENTED IIR NOTCH FILTER: P=3 


INPUT DATA FOR PROBLEM # 1 


PROBLEM # 1 INPUT DATA SOURCEFILE: diffeq.tst 
THE NUMBER OF INPUT DELAYS: L= 5 

THE NUMBER OF OUTPUT DELAYS: N= 5 
THE VALUE OF nstop IS: 400 
THE COEFFICIENTS b(0), b(1), ..., bCL) ARE: 


-990000E+00 -.329900E-02 .186124E-01 .937440E-02 =. 661007E+00 -.60470GRaae 


Mig COBPFICIENES saa), ., @(N) ARE: 


-OQQ000E+00 .O00000E+00 .0000O0E+00>.652382E1CO a3 ce59E+00 


OUTPLT DATA FOR PROBLEM # 1 


ns xCrs) y(ns) 

-5 . QOOCO0E+00 . OOOCOODE+00 
-4 . COOOOCE+00 . COOOD0E+00 
-3 . OOOOO0E+00 . OOOOO0E+00 
-2 . QOOOO0E+00 . OCOOOO0E+00 
-1 . QOOOOQ0E+00 . OOCOO0E+00 
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